o 

(N 

< 



< 



(N 
> 

(N 

O 
^, 

en 

o 

(N 



% 



INVARIANT MEASURES AND THE SOLITON 
RESOLUTION CONJECTURE 

SOURAV CHATTERJEE 



Abstract. The soliton resolution conjecture for the focusing nonhnear 
Schrodinger equation (NLS) is the vaguely worded claim that a global 
solution of the NLS, for generic initial data, will eventually resolve into 
a radiation component that disperses like a linear solution, plus a lo- 
calized component that behaves like a soliton or multi-soliton solution. 
Considered to be one of the fundamental open problems in the area 
of nonlinear dispersive equations, this conjecture has eluded a proof or 
even a precise formulation till date. 

This paper proves a "statistical version" of this conjecture at mass- 
subcritical nonlinearity, in the following sense. The uniform probability 
distribution on the set of all functions with a given mass and energy, if 
such a thing existed, would be a natural invariant measure for the NLS 
flow and would reflect the long-term behavior for "generic initial data" 
with that mass and energy. Unfortunately, such a probability measure 
does not exist. We circumvent this problem by constructing a sequence 
of discrete measures that, in principle, approximate this fictitious prob- 
ability distribution as the grid size goes to zero. We then show that a 
continuum limit of this sequence of probability measures does exist in a 
certain sense, and in agreement with the soliton resolution conjecture, 
the limit measure concentrates on the unique ground state soliton. This 
is proved by first proving a similar result for the discrete NLS. 

The above result, following in the footsteps of a program of study- 
ing the long-term behavior of nonlinear dispersive equations through 
their natural invariant measures initiated by Lebowitz, Rose and Speer, 
and carried forward by Bourgain, McKean, Tzvetkov, Oh and others, is 
proved using a combination of techniques from large deviations, PDE, 
harmonic analysis and bare hands probability theory. It is valid in any 
dimension. 
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1.1. Probabilistic motivation. Suppose that we are asked to choose a 
function / : M"^ — )• C uniformly at random from the set of all v : W^ — )• C 
satisfying AI{v) = m for some given constant m, where 

(1) M{v) := I \v{x)\'^dx. 

While this question does not make sense mathematically, the only reasonable 
answer that one can give is that / must be equal to zero almost everywhere. 
Paradoxically, this / does not satisfy M{f) = m. The paradox is resolved if 
we view this question as the limit of a sequence of discrete questions: First 
approximate M by a large box [—L, L] ; then discretize this box by splitting 
it as a union of many small cubes; finally, choose a function / : M — )• C 
uniformly from the set of all functions v that are piecewise constant in these 
small cubes and zero outside the box [—L, L] , and satisfy M{v) = m. This 
is a probabilistically sensible question; the resulting / approaches zero in 
the L°° norm as the box size goes to infinity. 

Now suppose that we add one more constraint, namely, that / should 
satisfy H{f) = E, where H is the functional 

(2) H{v):=- f \Vv{x)\'^dx ^— f \v{x)\'P+^dx, 

and p > 1 and E £ M. are given constants. The motivation for adding this 
second constraint comes from the study of microcanonical invariant mea- 
sures of nonlinear Schrodinger equations (more on this later) . One problem 
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that arises immediately is that if v satisfies M[v) = m and H{v) = E, so 
does the function u{x) := aov{x + xq) for any xq G M*^ and ao ^ S^, where 
S^ is the unit circle in C Thus, it is reasonable to first quotient the function 
space by the equivalence relation ~, where u ~ u means that u and v are 
related in the above manner. 

When p satisfies the "subcriticality" condition p < 1 + 4/d, standard re- 
sults from the theory of nonlinear Schrodinger equations imply that the set 
of functions v that minimize H{v) given M{v) = m form a unique equiva- 
lence class of the relation ~. This equivalence class is known as the "ground 
state soliton" of mass m. The main result of this manuscript (Theorem [T]) 
says that if we attempt to choose an equivalence class uniformly at random 
from all classes satisfying M(v) = m and H{v) = E, by first discretizing the 
problem and then passing to the continuum limit, then we end up choosing 
this ground state soliton. As before, there is no paradox in the fact that the 
ground state soliton may not satisfy the constraint H{v) = E. While the 
problem is quite simple for the single constraint M(v) = m, the addition of 
the second constraint H{v) = E somehow renders it unreasonably difficult; 
indeed, nearly the entirety of this long manuscript is devoted to the proof 
of Theorem [TJ 

The above result is a small step towards understanding uniform proba- 
bility distributions on manifolds in function spaces that are defined by a 
finite number of constraints. These distributions arise as "microcanonical" 
invariant measures for Hamiltonian flows on such manifolds. The conserved 
quantities for the flow give the constraints defining the manifold. A prelimi- 
nary attempt with a simpler problem was made in p^. All of this, and how 
it connects to the behavior of nonlinear Schrodinger flows and ideas from 
statistical physics, will be discussed in greater detail in the remainder of this 
section. 

1.2. The nonlinear Schrodinger equation. A complex-valued function 
u of two variables x and t, where x G M is the space variable and t G M is 
the time variable, is said to satisfy a d-dimensional nonlinear Schrodinger 
equation (NLS) if 

(3) idtu = —Au + k\u\p~^u, 

where A is the Laplacian operator in R , p > 1 is the nonlinearity parameter, 
and K is a parameter which is either -|-1 or —1. When k = 1, the equation 
is called "defocusing" , and when k = — 1 it is called "focusing" . 

The study of the NLS and other nonlinear dispersive equations is a large 
and growing area in the analysis of PDE, with numerous open questions and 
conjectures. For a very readable general introduction, see Tao [62]. For a 
more specialized account of the state of affairs in the study of NLS, see the 
lecture notes of Raphael [47j . The NLS arises in many areas of the pure and 
applied sciences, including Bose-Einstein condensation, Langmuir waves in 
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plasmas, nonlinear optics, and a number of other fields [2 H [22 | [30 } [7T | H} [23 | 
[7Ql[52]. 

The NLS is an infinite dimensional Hamiltonian flow, with Hamiltonian 
given by 

H{v) = l [ \Vv{x)\'^dx + ^^ [ \v{x)\P+Ux. 

(Note that in the focusing case {k = —1), this is just the function H defined 
in ([2]).) Consequently, if u is a solution to Q, then H{u{t,-)) is the same 
for all t. Since H{v) is commonly called the energy of v in the context of 
Hamiltonian flows, the previous sentence can be restated as: "The NLS flow 
conserves energy" . Another important conserved quantity is the mass M{v), 
defined in ([I]). 

A significant amount of information is known about the defocusing NLS; 
in particular, it is known that in many situations, solutions of the defocusing 
equation disperse like solutions of the linear Schrodinger equation (see [62| 
p. 154]). Here "dispersion" means that while M(n(t,-)) remains conserved, 
for every compact set K (^M. , 

lim / 1^(2;, t) I dx = 0. 

In the focusing case, however, dispersion may not occur. This is demon- 
strated quite simply by a special class of solutions called "solitons" or "stand- 
ing waves". These are solutions of the form u{x,t) = t;(x)e"^*, where cj is a 
positive constant and the function v is a solution of the soliton equation 

(4) -ojv = -Av -\v\P~^v. 

Often, the function v{x) is also called a soliton. Such functions are known 
to be smooth and exponentially rapidly decreasing (see e.g. [131 Section 
8.1]), and if one makes the further assumption that v is non- negative and 
spherically symmetric then there is a unique solution to ([!]) for each cj > 
[ini EOl [3]; we refer to this v as the "ground state". There also exist 
radial solutions which change sign, see [2j; such solutions are called "excited 
states" . 

The focusing equation is said to have mass-subcritical nonlinearity if the 
nonlinearity parameter p satisfies the subcriticality condition 

4 
1 <p < 1 + -. 
a 

Mass-subcritical nonlinearity has important consequences. For instance, if 
p < 1 + A/d, then it is easy to show that all solutions with initial data in H^ 
are global and bounded in H^ (see |47[ Section 1.1]). Another important 
feature of this regime is that for any m > 0, 

(5) S„,in(m):= inf H{v) G {-oo,0), 

V : M{v)=m 
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and the infimum is achieved at the ground state sohton with mass m. In 
fact, it is simple to prove by a scaUng argument that when p < 1 + 4/d, the 
function -Emin has the form 

(6) £^mm("^) = m°'E^in{l), 

where ii^min(l) £ (— co,0) and a > 1 is a constant that is exphcitly deter- 
mined by p and d (see [37., Section 1.4]). The infimum is achieved uniquely: 
any energy minimizing function v must be of the form 

v{x) =Qx(m){x-xo)e''^°, 

where xq S W^ and 70 G M, and Q\{rn) is the unique ground state soliton 
with mass m. The ground state soliton Q\Un) has the following explicit 
form: Let Q be the unique positive and radially symmetric solution of the 
equation 

(7) -Q = -AQ-|Q|P-iQ. 
For each A > 0, let 

(8) Qa(x) := A2/(p-i)Q(Ax). 

Then for each tti > 0, there is a unique \{m) > such that Qx{m.) is the 
ground state soliton of mass m. For all of the above claims about ground 
state solitons in the mass-subcritical regime, see \47\ Sections 1.2 and 1.3]. 
The uniqueness of the ground state is a particularly deep result. See also 
[62| Appendix B]. 

When p > 1 + 4/(i, much less is known; it is currently an area of active 
research (see [261 [2H] for recent developments and pointers to the literature) . 

Even in the mass-subcritical case, little is known about the long-term be- 
havior of solutions. One particularly important conjecture, sometimes called 
the "soliton resolution conjecture" (see Tao [62l p. 154]), claims (vaguely) 
that as i — )■ 00, the solution u{-,t) would look more and more like a soliton, 
or a union of a finite number of receding solitons. The claim may not hold 
for all initial conditions, but is expected to hold for "most" (i.e. generic) 
initial data. In the critical and supercritical regimes, the conjecture is still 
supposed to be true, but with the additional imposition that the solution 
does not blow up. The conjecture is based mainly on numerical simulations, 
although there has been a limited amount of progress towards a proof (see 
[401 [56| \6T\ l63| [64] and references therein) . The only case where one can give 
a heuristic treatment is when d = 1 and p = 3, where the NLS is completely 
integrable (see [55l|4Tl[72]). The soliton resolution conjecture has been inves- 
tigated for other dispersive systems, with partial results [34 l [53 t [54iri9| l20] . 
For recent progress on the soliton resolution conjecture for the energy-critical 
wave equation and a far more extensive survey of the literature around the 
conjecture, see [18] . 
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1.3. Invariant measures for the NLS. One approach to understanding 
the long-term behavior of global solutions is through the study of invari- 
ant Gibbs measures. Roughly, the idea is as follows. Since the NLS is a 
Hamiltonian flow, one might expect by Liouville's theorem that Lebesgue 
measure on the space of all functions of suitable regularity, if such a thing 
existed, would be an invariant measure for the flow (see, e.g., p?, P- 68] for 
a statement of Liouville's theorem in the finite dimensional setting). Since 
the flow preserves energy, this would imply that Gibbs measures that have 
density proportional to 

(9) exp{-Kl3H{v)) 

with respect to this fictitious Lebesgue measure (where /3 is arbitrary) would 
also be invariant for the flow. One way to make this rigorous is to first 
restrict the system to the unit torus T'^ and then consider Gibbs measures 
that have density proportional to 

(10) exp(K/3 / \v{x)\P^^dx 

with respect to the free-field Gaussian measure (see [32]) on the appropriate 
space of distributions on T'^. This is the pioneering idea of Lebowitz, Rose 
and Speer [33]. For such a thing to make sense in d > 2, one has to interpret 
the integral in the Wick-ordered sense. 

These Gibbs measures exist for the defocusing case (k = 1) for all p in 
d = 1 (without Wick ordering) and for p < 5 in d = 2, and p < 3 in 
d = 3 [23]. Furthermore, despite the fact that this measure is supported on 
rough functions, Bourgain showed that it is invariant under the dynamics 
given by ^ for d < 2 [6]. This means that the dynamics can be defined 
(after Wick ordering modification in d = 2) on a set of full measure with 
respect to this Gibbs measure. 

The focusing case {k = —1) is more delicate. Since H is unbounded 
from below, it is obvious that the Gibbs measure cannot exist without some 
restrictions on its domain. It was shown in [33] that in d = 1, the Gibbs 
measure exists for p = 3 when restricted to L^ balls and that it exists for p = 
5 with the additional condition of small /3. The development was continued 
by Bourgain [5] , McKean [35] , McKean and Vaninsky in [Ml EZj [38] and 
Zhidkov [73] . In d = 2, Jaffe showed that the measure exists for p = 2 for 
real u when restricted to L^ balls and after Wick ordering (see [32]); while 
Brydges and Slade [9] showed that this does not work when p = 3. 

Invariant measures coupled with Bourgain's development [HIHllTllH] of the 
so-called X*'* spaces ("Bourgain spaces") for constructing global solutions 
has led to important developments in this field. Recently, striking advances 
have been made by Tzvetkov and coauthors [Ml [13 [671 [HI [El EH IM] and 
Oh and coauthors [iTl [Ml |43l El [Ml [IS] and others (e.g. [39]) who use 
invariant measures and Bourgain's method to construct global solutions of 
the NLS and other nonlinear dispersive equations with random initial data. 
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Qualitative features of the infinite volume limit of Gibbs measures were 
studied by Brydges and Slade [9] in d = 2 and Rider [48l [49] in d = 1. In- 
variant Gibbs measures for the cubic discrete nonlinear Schrodinger equation 
(DNLS) in d > 3 were studied in [TS] . 

Another idea that is gaining traction in the physics circles in recent years 
is that of considering micro canonical ensembles (see e.g. [511 [52] and refer- 
ences therein). The general idea ~ which has already been discussed at the 
beginning of this section - is to consider an abstract manifold of functions 
satisfying certain constraints (usually two) and then trying to understand 
the characteristics of a function picked uniformly at random from this man- 
ifold. Often, the physicists alternately characterize the uniform distribution 
as the "maximum entropy" distribution. For example, in the context of the 
NLS, one looks at the "uniform distribution" on the space of all functions 
with a given mass and energy. The relevance of this to the long term be- 
havior of NLS flows is heuristically justified through Liouville's theorem; we 
have more on this in the next section. It is in an attempt to understand these 
physical heuristics that I got interested in this line of research (and I thank 
Persi Diaconis - who heard about it from Julien Barre - for communicating 
these problems to me a few years ago) . In an early paper [T3] that I wrote on 
this topic, I tried to understand the behavior of functions chosen uniformly 
from all functions satisfying f\u{x)\'^dx = m and j \u{x)Y'~^^ dx = —E, 
completely ignoring the gradient term in the Hamiltonian. Already in this 
simplified situation one can prove interesting phase transitions and localiza- 
tion phenomena. In a later paper with Kay Kirkpatrick [15], the gradient 
term was added to the analysis, but the nonlinearity parameter was taken to 
be so large that the gradient term became practically unimportant. In both 
[H] and [15], the settings were discrete and too crude to allow passage to 
a continuum limit. The purpose of the current manuscript is to undertake 
the more serious task of analyzing regimes where the gradient term actually 
matters, and a continuum limit can be taken. 

2. The main result 

Assume that k = —1 for the rest of this manuscript. Given tti > and 
E > £'inin(7n) (where £'niin("i) is the minimum energy for mass m, as defined 
in ([SD), let 

(11) S{E, m) := {v £ F^(M'^) : M{v) = m and H{v) = E} 

be the set of all H^ functions of mass m and energy E. Since the NLS 
flow ([3|) preserves mass and energy, the same heuristic via Liouville's theo- 
rem that led to ([9|) would imply that a "uniform distribution" on S{E,m), 
if such a thing existed, would be an invariant measure for the flow. In 
physics parlance, these measures would be the "Microcanonical Ensembles" 
corresponding to the "Grand Canonical Ensembles" given by ([9|). 

If the soliton resolution conjecture is indeed true, and an invariant mea- 
sure like the microcanonical ensemble suggested above indeed exists and 
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describes the long-term behavior of the typical NLS flow with a given mass 
and energy, then it should put all its mass on soliton or multi-soliton func- 
tions. This may seem like a contradiction since such functions may not have 
energy E as required for membership in S{E, m). However there is no actual 
contradiction since S{E, m) is not compact under any reasonable metric. 

Our goal is to go ahead and try to give a meaning to the abstract nonsense 
outlined above. To give a meaning to the notion of a uniform probability 
distribution on the set of all functions with a given mass and energy, we 
restrict ourselves first to a finite region of space, and then to a discretization 
of it. Instead of M"^, therefore, our space would be the discrete grid 

K = {0, 1, . . . , n - 1}^ = {Z/nZ)"^. 

We imagine this set embedded in M as hVn, where /i > is a parameter 
representing the grid size. Note that hVn is a discrete approximation of the 
box [0,nh] . We would eventually want to send nh to oo. 

The mass and energy of a function v : Vn ^ C at grid size h and box size 
n are defined in analogy with ([T]) and ([2]) as 

(12) Mi^M ■■= h" V Hx)\\ 



xeVn 



and 



/i'^ 



(13) HhAv)--='-^ E 



v{x) - v{y) 



h 



2 



^\ E 1^(^)1""''' 



V + - 



where x ~ y means that x and y are neighbors in Vn- For simplicity, we 
endow Vn with the graph structure of a discrete torus, i.e. identifying n 
with 0. Let C^" denote the set of all functions from Vn into C Take any 
e > 0, rn, > and £^ G M, and define the set 



(14) S,^h,n{E, m) := {v G C^" : \Mh,n{v) - m\ < e, \Hh,n{v) - E\ < e}. 

Clearly, S^^h^n{E^m) is a finite volume subset of the finite dimensional space 
C^". This set is a "manageable" version of the set S{E,m) defined in pi|) . 
Indeed, as e — )• 0, /i — )• and nh ^ oo, the set S^^h,n{E, m) may be imagined 
as tending to the limit set S{E,'m). 

We have chosen e > to ensure that the volume is nonzero whenever the 
set is non-empty. In this situation, the uniform probability distribution on 
this set is well-defined. This uniform probability distribution, besides being 
an approximation to our abstract object of interest, has also a concrete 
interpretation as a natural invariant measure for an appropriate discrete 
NLS evolution on Vn, to be discussed in the next section. 

Fix E and m such that E > E^i^{m). Given e, h and n, let ^e,h,n be 
the uniform probability distribution on Se,h,n{E, m). Let fe,h,n be a random 
function on Vn with law He,h,n- Our main result, stated below, is that when 
the nonlinearity is mass-subcritical, the random function fe,h,n converges in 
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a certain sense to the unique ground state soliton Qx{ni) of mass m defined 
in Section [H as {e,h,nh) goes to (0,0, oo) in a certain manner. 

To define the notion of convergence, we first need a way of comparing 
functions on Vn with functions on Z and M . Given v : Vn ^ C, first 
define its extension v'^ to Z by simply defining v^ as equal to t; on V^ and 
zero outside. Next, given a function w : Z^ — t- C, define its "continuum 
image at grid size h" as the function w -.W^ ^ C, defined as follows. Given 
y = {yi, . . . , yd) G M'^, let x = (xi, . . . , Xd) be the unique point in Z'^ such 
that for each i, 

Xi < yi/h < Xi + l, 
and let w{y) := w{x). (In other words, Xi = [yi/h].) Lastly, given v : Vn ^ 
C, define its continuum image i) : R"' — )• C at grid size h as the function v'^, 
that is, the continuum image of the extended function v^. 

For each q G [l,oo] define a pseudometric L' on the set of measurable 
complex- valued functions on R as 

(15) L'^{u,v) := inf \\u{-) - aov{- + Xo)\\q, 



where S^ is the unit circle in the complex plane and \\-\\q denotes the usual L"? 
norm of a complex- valued function on R"^ with respect to Lebesgue measure. 
This is a pseudometric since L'^{u,v) may be zero even if u and v are not 
equal, but v is of the form v{x) = aou{x + xo) for some xq G M"^ and oq G S^. 
It is necessary to work with pseudometrics since the law of fe,h,n is invariant 
under translations and multiplication by scalars of unit modulus. 

Theorem 1. Suppose that 1 < p < 1 + 4/d. Fix E and m such that 
E > £'inin(m), where E'min is defined in ([5]). Let fe,h,n be a uniform random 
choice from the set S^^h,niE,m) defined in ([HI), and let fe,h,n denote its 
continuum image at grid size h, as defined above. Let L'^ he the pseudometric 
defined above. Then for any 6 > Q and any q G (2,cxd], 

lim limsuplimsupP(L«(/,^;,_„,(5;,(.^)) > 6) = 0, 

h-^O e_!.o n-s>oo 

where Qxun) is the unique ground soliton of mass m defined in Section [H 
Moreover, there is a sequence (efc,/ifc,nfc) with e^ — )■ 0, /i^ — )■ and n^hk -^ 
oo as k ^ oo, such that for any fixed 6 > and q G (2, cxd], 

limP(L«(4,;,,,„,,Q;,(„))><5)=0. 

Note that the energy of fe.h,n is converging to E, whereas the energy 
of Q\(m) is Ejnini'm-), which is strictly less than E. There is no contradic- 
tion here, since the metrics of convergence are not strong enough to ensure 
convergence of the Hamiltonian. Nor should they be, actually, since the 
difference E — E^i^{m) denotes the amount of energy that has "escaped" to 
infinity when the NLS has flowed for a long (infinite) time. 

Let me emphasize here that Theorem [T] does not say that a certain Gibbs 
measure concentrates on its lowest energy state in the infinite volume limit. 
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In fact, the soliton Q\(^rn) is not in the support of the measure iJ,e,h,n at all, 
even in the limit. What Theorem [T] says is much more subtle: In the infinite 
volume continuum limit, a typical function with mass m and energy E de- 
composes into an "invisible" or "radiating" part that is small in L°° norm 
but contains a significant amount of energy due to microscopic fluctuations, 
and a "visible part" that is close to the soliton Q\Un) in the L°° distance. 

It may seem strange that while a certain amount of energy escapes to in- 
flnity, there is no escape of mass. Again there is no contradiction, since func- 
tions of arbitrarily small L^ mass on M can hold arbitrarily large amounts 
of energy by being very wiggly. 

Lastly, one may question the implications of Theorem [1] for the dynamics 
of the NLS. Indeed, a probability measure that concentrates on the time 
oscillates of a single soliton is also an invariant measure for the NLS and 
carries no special information. The difference between such a special in- 
variant measure and the Microcanonical Ensemble is that support of the 
Microcanonical Ensemble contains almost every function of mass m and en- 
ergy E. In other words, almost every ergodic component of functions with 
mass m and energy E has the property that for initial data in that ergodic 
component, the flow converges to the soliton Qum) as time tends to infinity. 

Of course, none of this makes rigorous sense since we are in a highly 
non-compact setting. In fact, it is clear that Theorem [1] does not model 
the full dynamics of NLS from the following observation: There are stable 
multisoliton solutions consisting of two or more receding solitons which do 
not generally collapse into a single ground state. This is possibly because the 
effect of recession "outruns" the thermodynamic convergence to equilibrium. 

Nevertheless, the argument leading to Theorem [1] does give a proof of a 
statistical version of the conjecture in the discrete setting when we take e 
very close to and n very large, but keep h fixed. In this situation there 
is enough compactness to ensure that receding solitons cannot outrun the 
convergence to equilibrium. This is the topic of the next section. 

3. The soliton resolution conjecture for the discrete NLS 

Equation ([3]) can be discretized to give the so-called 'discrete nonlinear 
Schrodinger equation' (DNLS). For the DNLS, only space (and not time) is 
discretized: we substitute W^ by the discrete grid /iZ"', where h is the lattice 
spacing. The limit /i — t- is called the 'continuum limit' in the literature on 
the DNLS. 

The discrete nearest-neighbor Laplacian on Z"' with grid size h is de- 
fined as 

(16) Av{x):=^ Yl iHy)-v{x)), 

y ■■ y~x 

where y ~ x denotes the sum over all neighbors of x and v is any complex- 
valued function on Z'^. Note that the scaling by h"^ is meant to ensure 
that if V is the restriction of some smooth function v on the grid /iZ (i.e. 
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v{x) = v{hx)), then the discrete Laplacian converges to Av as the grid size 
goes to zero. 

In analogy with equation ([3]), the discrete focusing NLS on hZ, with 
nonhnearity parameter p is a family of coupled ODEs 

(17) i^ = -Au-\u\P-\ 
^ ^ dt II' 

where A is the discrete Laplacian for grid size h defined in (J16p and u{x, t) 
is a function on Z x M. The mass Mfi{v) and the energy Hh{v) of a function 
u : Z'^ — )• C at grid size h may be defined just as in ([12]) and (fT3]) . but after 
replacing Vn by Z . As in the continuous case, the DNLS is a Hamiltonian 
system and the mass and energy are conserved under the DNLS flow. 

The DNLS has been studied widely by physicists, but not so much by 
mathematicians, particularly in dimensions higher than 1. For a recent 
survey of the mix of rigorous and non-rigorous results that exist in the 
literature, see [27j. We will now outline some of the results about the DNLS 
that are proved in this manuscript en route to the proof of Theorem [TJ 

Given /i > and m, > 0, let E^axim,h) and E^\n{m,h) denote the 
supremum and infimum of possible energies of functions with mass m at 
grid size h. 

Theorem 2. Suppose that 1 < p < 1 + 4/d. With the above definitions, 
for any h > and m > we have E^n;,^{m,h) = 2dm/h'^ and — cxd < 
Erain{fn,h) < 0. Moreover, the function -Emin satisfies for all positive m 
and m' the strict subadditive inequality E-^i^{m + m',h) < E^i^{m,h) + 
Emm{rn',h). Lastly, linih^Q Ejamiin, h) = Ejainiin) and the convergence is 
uniform over compact subsets of {0,oo). 

The first couple of assertions of Theorem [2] are proved in Section [HI The 
subadditive inequality is proved in Section [T5j The convergence argument is 
more complicated. It follows from Corollary [69] in Section [18] The conver- 
gence is based on the convergence of discrete solitons to continuum solitons 
(see Theorem [5] below) . These scattered results are gathered into a formal 
proof of Theorem [2] in Section [20] 

As in Section [2] we define a set of pseudometrics on the space of all 
complex- valued functions on Z*^. For any q E [l,oo], let L^ be the pseudo- 
metric on C^ defined as 

(18) L'i{u,v) := inf \\u{-) - aQv{- + XQ)\\q, 



where || • ||g is the usual L^ norm on C . 

Let S{m,h) be the set of all functions / with Mh{f) = m and Hh{f) = 
EminiiT^, h). The set S{m, h) will be called the set of discrete ground state 
solitons with mass m at grid size h. Note that, as in the continuum case, 
a simple Euler-Lagrange argument shows that any discrete ground state 
soliton must necessarily satisfy the discrete soliton equation 

—iov = —Av — \vF~ v 
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for some w > 0. 

Unlike the continuum case, the discrete ground state sohton for a given 
mass may not be unique. However, we do know from the following theorem 
that S{m, h) is non-empty and compact in the L'^-topologies. Not only that, 
the set S{m, h) also has an analog of the so-called "orbital stability" property 
(see [471 Section 1.3]) of continuum solitons: any function that has near- 
minimal energy must be nearly a soliton. The subadditive inequality from 
Theorem [21 together with the classical concentration-compactness technique 
[47[ Section 1.4], is the key to the proof of this result. (Note that the 
subadditive inequality is trivial in the continuous case by the formula dG]).) 

Theorem 3. Suppose that 1 < p < 1 + 4/d. Let S{m, h) he the set of 
ground state solitons of mass m at grid size h, as defined above. Then for 
any m > and h > 0, S{m,h) is non-empty. Moreover, for any sequence 
of functions fk such that Mh{fk) — >• m and Hh{fk) — )• £'jnin(m, h), there is a 
sub-sequence fk and some f G S{m,h) such that fk converges to f in the 
L^ pseudometric for any q G [2,oo]. 

The main argument for the proof of Theorem [3l using a discretization 
of the concentration-compactness method, is presented in Section [T5j The 
proof is completed in Section [20l 

What happens to S{m, h) as h tends to zero? The next theorem answers 
this question. As /i — >• 0, the S{m, h) shrinks to a single point, namely, 
the unique continuum ground state soliton Q^m) defined in Section [H For 
related results on continuum limits of the discrete NLS in one dimension, 
see [29]. 



Theorem 4. Suppose that 1 < p < 1 + 4/d. Let ruk be a positive se- 
quence converging to some m > 0. Let hk be a positive sequence tending 
to zero. For each k, let fk be an element of S{mk,hk). Let fk be the con- 
tinuum image of fk at grid size h, as defined in Section [2j Let L'^ be the 
pseudometric on L'^iW^) defined in Section^ Then for any q G [2,cx)], 

lilTlfc^oo L'i{fk, Q\{m)) = 0. 

The above theorem is proved by showing that for small /i, discrete ground 
state solitons can be approximated by smooth functions - and then applying 
the orbital stability of continuum solitons [471 Section 1.3]. To show that 
discrete ground state solitons can be approximated by smooth functions, 
one has to prove regularity estimates that do not blow up as /i — )• 0. To 
achieve this, the route taken in this paper is to translate the proof of reg- 
ularity of continuum solitons (as sketched in |62[ Proposition B.7]) to the 
discrete setting and obtain "/i-free" estimates. This translation necessitates 
the development of a slew of discrete harmonic analytic results, including 
fine properties of discrete Green's functions, discrete Littlewood-Paley de- 
compositions, and a discrete Hardy-Littlewood-Sobolev inequality of frac- 
tional integration. The harmonic analytic tools are developed in Section [T6| 
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and the regularity of discrete solitons is worked out in Section [T71 The 
convergence argument is presented in Section [HI 

Next, recall the random function fe,h,n from Section [2j Observe that by 
Liouville's theorem, the law of fe,h,n is an invariant measure for the DNLS 
flow (fTTj) at grid size h. In Theorem [H we saw what happens to fe,h,n as 
(e, h, nh) — )• (0, 0, oo) in a certain manner. On the way to proving Theorem[Tl 
we first investigate what happens to fe,h,n as (e,n) — )• (0, oo), fixing h > 0. 
What happens is the following: with high probability, fe^h,n is close to a 
discrete ground state soliton of mass m' , where m' S [0,m] is determined in 
a complicated manner hy E, m and h. The following theorem makes this 
precise. 

Theorem 5. Suppose that 1 < p < l+A/d. Take any E €M., m > and h > 
such that E^[^{m,h) < E < E^n_^{m,h). Let fe^h,n be a uniform random 
choice from the set S^^h^n{E, m) defined in (J14p . and extend its domain to U^ 
hy defining it to he zero outside Vn- If E < i£'max("i-, /i), then there exists 
a compact set K C [0, m] such that for any 5 > Q and any q S (2, oo], 



lim lim P inf inf L'^ifehn, v)>6] =0. 

e-s-0n-5>oo \m'£KveS{m',h) ' ' / 

Furthermore, the set K can he described as follows. It is the set of all 
m! G [0, m] that maximize 

log(m - m ) - *d , 

\ m — m J 

where ^d : M — t- [0, oo] is the function 

(19) ^d(a) = sup / logfl — 7H y^ siT?{-KXi)\dxi- ■ ■ dxd 

0<7<1 J[o,i]<* V ^ i^l / 

for a G (0, 2d), ^d(a) = ^d(4d - a) for a G (2d,4d), ^di'^d) = 0, and 
'^d{ot) = oo for a > Ad and a < 0. Lastly, if E > ^E^^.^^{m,h), then for 
any (5 > 0, 

lim lim P(||/,,mIIoo ><5) = 0. 

In a few words, the above theorem says the following: "For the DNLS on 
a large torus, a typical function with a given mass and energy is close in the 
L°° distance to a soliton with a (possibly) different mass and energy." 

Theorem [5] proves the following statistical version of the soliton resolution 
conjecture for the DNLS restricted to a large torus. The uniform distribu- 
tion iJLe,h,n on the set Se,h,niE, m) is an invariant measure for the DNLS flow 
on the torus Vn but is not ergodic. Elementary ergodic theory tells us that 
it decomposes as a convex combination of ergodic components (see e.g. [69^ 
Section 6.2]). The ergodic decomposition of fie,h,n yields a natural probabil- 
ity measure Te^h,n on the space of ergodic measures on Se,h,n{E,m) for the 
DNLS flow (that chooses ergodic components with "probability proportional 
to size"). 
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The long-term behavior of the flow is dictated by the ergodic component 
to which the initial data belongs. Theorem [5] implies that if T^^h,n is the 
relevant probability measure on the space of ergodic components, then for 
the vast majority of components, the long-term behavior indeed matches 
the prediction of the soliton resolution conjecture: some amount mass and 
energy flow out "to infinity" in the form of a function of small L°° norm, 
and whatever remains is close to a soliton. 

The proof of Theorem [5] is divided into several components. The first 
ingredient, proved in Section O is a large deviation principle for gradients of 
random functions. The variational problem related to this large deviation 
principle is dealt with in Section [TOl The proof of the main identity in 
Theorem [5] is separated into two pieces: first an upper bound, and then a 
matching lower bound. The upper bound in the case E < ^-Emaxl^i) h) is 
proved in Section [TTJ The matching lower bound requires us to first prove 
the exponential decay of discrete solitons. This is done in Section [12l mostly 
along the lines of the proof of exponential decay of continuum solitons |62l 
Proposition B.7], but with the crucial difference that we now have to deal 
with discrete Green's functions. Using the information from Section \T2\ the 
lower bound is proved in Section [T3j Finally, the case E > 2-E'max("i-) ^) is 
handled in Section [TJl Everything is formally put together to complete the 
proof of Theorem [5] in Section [JH 

Last of all, let us indicate how all this leads to the proof of Theorem [TJ 
This is quite easy, given Theorems [5] and HI We just take the limit /i — )• 
in the explicit formula given in Theorem [SJ The main step is to prove 
that the set K in Theorem [5] shrinks to the singleton set {0} as h goes 
to 0. This is done in Section [191 The convergence of discrete solitons to 
continuum solitons as given by Theorem H] finishes the proof. This argument 
is formalized in Section! 



4. Main ideas in the proof 

Let / : M — 7- C be a function "uniformly chosen" from the set of functions 
satisfying M(f) = m and H{f) = E, whatever that means. We need to show 
that for any set A of functions that do not contain the ground state soliton, 
the chance oi f € A is zero. 

Take any 6 > and let Vs := {x : \f{x)\ < 6}. Then 

\f{x)\P+^dx < 6P-^ [ \f{x)\^dx < 6P-^m. 
Vs Jvs 

Decompose f as u + v, where u = flvg and v = fl^d\y^. The above 
inequality shows that when 5 is close to zero. 



H{u) « - / \Vu{x)\'^dx. 



INVARIANT MEASURES AND THE SOLITON RESOLUTION CONJECTURE 15 

On the other hand 

Vo1(mV5)<4 / l/(x)p(ix<^. 



Let us refer to v and u as the "visible" and "invisible" parts of /. The last 
two inequalities show that: 

• The visible part is supported on a finite volume set, whose size is 
controlled by 5. 

• The energy of the invisible part is essentially the same as the L^ 
norm squared of its gradient, times 1/2. 

The game now is to compute P(/ G A) by controlling the visible and invis- 
ible parts separately. The visible part, being supported on a "small" set, 
can be analyzed directly. For the invisible part, one has to develop joint 
large deviations for the mass and the gradient, since the nonlinear term is 
negligible in the invisible part. Solving the variational problem related to 
the large deviation question, one arrives at the conclusion that the visible 
part must be close to the ground state soliton with high probability. 
The main steps in the above program are the following. 

(1) Develop large deviation estimates for the invisible part in the finite 
volume discrete case. This is done in Sections El [71 [9l [TTl [121 [El 

and[Il 

(2) Analyze the variational problem related to this large deviation ques- 
tion, and thereby show that with high probability, the visible part 
has the minimum possible energy for its mass. This is done in Sec- 
tions [HI [lO] and [El 

(3) Pass to the infinite volume limit (keeping the grid size fixed) using a 
discretization of the classical concentration-compactness argument, 
and show convergence to discrete solitons. This is done in Sections 
[El and [201 

(4) Develop discrete analogs of harmonic analytic tools (Littlewood- 
Paley decompositions, Hardy-Littlewood-Sobolev inequality of frac- 
tional integration, Gagliardo-Nirenberg inequality, discrete Green's 
function estimates, etc.) to prove smoothness estimates for discrete 
solitons that do not blow up as the grid size — )• 0. This is done in 
Sections [El and [El 

(5) Use these smoothness estimates, together with the orbital stability 
of the ground state soliton, to prove convergence of discrete solitons 
to continuum solitons. This is done in Section [181 

5. Summary of notation 

In this section we summarize the notation that will be used repeatedly 
in this manuscript. Some of it has already been introduced, and some will 
be defined in later sections. The summary in this section is for the reader's 
convenience. 
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5.1. Spaces and norms. For a typical element x E M , we denote the 
ith coordinate of x by xi. The usual Euclidean norm of a vector x S M is 
denoted by |x|, while the l^ norm of x is denoted by |x|i. The same notation 
is used for norms of vectors in Z . 

The L'^ norms for functions on M and Z are defined as usual, and they 
induce the pseudometrics L'^ defined in (|15p and (llSp . The L'^ norm of a 
function v, whether on M or on Z , is denoted by \\v\\q. 

Sometimes, we use a slightly different L^ norm for functions on Z"' by 
combining the usual L'^ norm with the grid size h to get an L^' -norm: 

||,,|| /)'^/'3|U;|| 

II 'J\\q,h ■ — "' II ''llg- 

These norms will be used heavily in Sections [16] and [T71 

5.2. The discrete torus. Assuming that the dimension d is fixed, the torus 
Vn in Z"' is the set {0, 1, . . . , n — l}*^. We say two elements x and y in Vn 
are neighbors, and write x ^ y, \i\x — y\ = 1, where the difference x — y \s 
computed by subtraction modulo n in each coordinate. 

Sometimes, we use dVn to denote the boundary of the torus when con- 
sidered as a subset of Z'^ (without the toric graph structure) . In general dU 
denotes the boundary of a set C/ C Z"^ or [/ C V^. That is, dU is the set 
of points in U that are adjacent to some point outside U . Note that this 
boundary is different if U is considered as a subset of the torus Vn rather 
than as a subset of Z . Similarly, V^ denotes the set Ij \U when f7 is a 
subset of Z*^, whereas U'^ denotes Vn\U when U is considered as a subset of 
the torus Vn- 

5.3. Mass and energy. The mass of a function t; : M'^ — )■ C is defined as 

M{v) = I \v{x)\'^dx, 



and its energy, in the context of the NLS equation ([3]) with k = —1, is 
defined as 

H{v) = - I \Vv{x)\'^dx ^— I \v{x)\P+^dx. 

When u is a function on Z'^, its mass "at grid size /i" is defined as 



Similarly, the energy at grid size h is defined as 

(20) ^,.(.,=^5: ^M_^^^^K,,l 






IP+1 



For a function v on the torus Vn, the mass Mh,n{v) and the energy Hh^n{v) 
are defined exactly as in the last two displays, except that the sums are now 
over elements of Vn instead of Z . 
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For a function f : Z — )• C, the energy Hh{v) can be decomposed into the 
gradient component Gh{v) and the nonhnear component Nh{v), defined as 
the first and second terms on the right-hand side in (|20p . so that Hh{v) = 
Gh{v) — Nh{v). Similarly, we define Gh,n{v) and Nh^niv)- 

Given a subset C/ C V^, we define 

h'^ -^-^ v{x) — v{y) 



Gh,n{v,U) = - ^ 

xr^y 

and Hh^n{v-,U), Nh^n{v-,U) and Mh^n{v,U) are defined similarly. 

5.4. Maximum and minimum energies. Suppose that p and d are given. 
For each m > 0, we define E^^^i^ and E^^i^{m) to be the supremum and 
infimum of the set of all possible energies for functions with a given mass on 
W^, where mass and energy are defined as above. It is not difficult to verify 
that £'max(fw) = CO, and it is known that E^[^{m) has the form ([6]) when 
p<l + A/d. 

In the discrete case on Z'^, we define Eyna.x{m,h) and £^min(m,/i) to be 
supremum and infimum of the set of all possible energies for functions with 
mass m, where both mass and energy are computed at grid size h. We prove 
a number of things about these quantities in Section [H 

Finally, for functions on the torus Vn, we similarly define E^ia^{m,h,n) 
and -Emin("i5 h,n). 

Two related functions E^ and E^ are defined Section [TOl 

The set of functions on U^ with mass m at grid size h that minimize 
energy is denoted by S{m, h). 

5.5. Notation related to the variational problem. The function ^^ 
defined in (I19p is of fundamental importance in this manuscript. Two 
other functions, and 0, and two related subsets of M^, TZ{E,m,h) and 
M. {E, m, h), are defined in the beginning of Section [lOl All of these are used 
repeatedly in the manuscript. 

5.6. Continuum image of a function. Given a function u : Z'^ — )• C, its 
"continuum image at grid size /i" is defined in Section [21 but let us repeat the 
definition here. The continuum image at grid size h is a function {i : M'^ — )• C 
defined as follows. Given y = (yi, . . . , yd) E M*^, let x = {xi, . . . , Xd) be the 
unique point in Z'^ such that for each i, 

Xi < y-i/h < Xj + 1, 

and let v{y) := v{x). When u is a function on Vn, we define the continuum 
image by first defining the function to be zero on Z \Vn, and then defining 
the continuum image as above. 
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6. Comparison with a Gaussian function 

Let 4> = {(f){x))x^v„ be a collection of i.i.d. complex Gaussian random 
variables, with E((/)(a;)) = E(0(x)^) = and 

Recall the function fe^h,n and the set S^^h,n{E,m) defined in Section [2j The 
following basic lemma connects the properties of (j) with that of fe,h,n- 

Lemma 6. Take any m > and E &M. Let f = fe^h,n ctncl S = S^^h,niE, m) 
for simplicity. Then for any A (^ S, 



Proof. For any measurable A C C ", 

Now, \Mh^niv) — m\ < € when v £ S. Therefore, 

P(0 eS)= ^^^^ \, dv 

< (n'^/i'^7r-^)"'e-"'("'-^)Vol(S). 
Similarly, 

/■ „-n'^Mh,„{v) 

> (n'^/i'^7r-^)"'e-"'("^+^)Vol(A). 
This completes the proof. D 

7. Diagonalizing the Laplacian 

Let r = {T{x,y))x,yeV„ be the matrix defined as 

2d if X = y, 
^{x,y) = {-I if X ~y, 

in all other cases. 

The matrix F may be viewed as an operator acting on C " in the natural 
sense. The action of F on a function / : V^ — )• C will be denoted by Tf. 
Notice that F = — /i^A, where A is the discrete Laplacian on the torus Vn 
with grid size h, as defined in (|16p . and that F is a real symmetric matrix 
of order n . In this section, we will write down the spectral decomposition 
ofF. 
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For two functions u,v ^ C^", let {u,v) be the standard inner product, 

{u,v) := J2 u{x)v{x), 

X&Vn 



where v{x) is the complex conjugate of v{x). Notice that for any v, 

{v,Tv) = {Tv,v) = ^ \v{x) - v{y)\^ . 



x.yeVn 



We use the notation Xi to denote the ith. coordinate of a vector a; G M . 
Lemma 7. For each y = (yi, . . . , y^j € Vn, let py be the function 

Then the functions iPy)yeVn form a complete orthonormal system of eigen- 
functions of T, and the eigenvalue corresponding to py is 

d 
Xy := 4^sin^(7ryi/n). 
1=1 
Proof. To prove orthogonality, first notice that for any k ^ Z,, r := e^^Trfc/n 
is an nth root of unity, and hence 

n— 1 n— 1 < 

ri = \ 

\n if fc = 0. 



y^ ^i2nkj/n ^ y^ ^j _ J if fc / 0, 



j=0 j=0 

Thus, for any y,y' £Vn, 

n— 1 

{py,Py>)=n~'^ Y^ e'2^((?/i-3/i)^i+-+(?/d-?/:,)^d)/n 



xi,...,xa=0 








d n-1 








-".... > : ^"' 


^{y,- 


-y> 


^/n 


j=l xi=0 








'o ify/y', 








1 if y = y'. 









To show that py is an eigenfunction of F with eigenvalue Ay, note that for 
any x € Vn, 



^Pyix) = Y^ (Pyix) - Py{z)) 

z : z^x 
d 



i=l 

d d 

= -py{x) ^(e-?'^/" - e--^»/")2 = 4p,(x) J]sin2(7ry,/n). 

j=l i=l 

This completes the proof of the lemma. D 
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Let R G C "^ " be the matrix whose yth column is py, for each y ^ V . 
Note that i? is a unitary matrix. Let A G c^n.xy„ ^^^ ^^^ diagonal matrix 
whose yth diagonal element is \y. Then 

r = RKR*, 

where R* is the adjoint of R. This is the spectral decomposition of F. 

8. The set of possible energies for a given mass 

Let £'niin(w,, /i, n) and E„ii^x(m, h, n) be the minimum and maximum pos- 
sible energies at grid size /i of a function / : V^ — )■ C with mass ra. Since the 
map / I—)- Hh,n{f) is continuous and {/ : Mh,n{f) = f^} is a compact con- 
nected subset of C^", therefore for every E G [E^\^{m, h, n),Eyna.x{'>n, h, n)], 
there exists some fonVn with mass m and energy E at grid size h. Recall 
that £^min(w-5 h) and Ema.x{rn, h) are the infimum and supremum of the set 
of possible energies of functions with mass m on the full lattice l/, at grid 
size h. 

Lemma 8. For any m > 0, and any E G {E^i^{m,h),Ei^ax{m,h)), there 
is a function / : Z'^ — )■ C such that Mh{f) = m and H^^f) = E. Moreover, 
for any m> 0, any function on TL'^ with mass m and energy E (at grid size 
h) satisfies 

\E\ < C{p,d,h){m + m^P+^^/^). 

Proof. If Mh{f) = m, then for any x, |/(x)p < mh^'^. Therefore, 

i/(x)r+^ < {mh-'')(p~'y^\f{x)\^. 

Thus, 



2dh ^m + 



f^-d(p-l)/2^(p+l)/2 



P+l 

This proves the inequality. 

Next, take any two functions f,g € C with Mh{f) = Mh{g) = m > 0. 
If / = (7 or / = —g, then Hh{f) = Hh{g). Suppose that / / ^g. For 
each e G [0, 1], let /, := ^/ + (1 - e)g. Since M^if) = Mh{g) and / / -g, 
it is easy to see that fe is not the zero function for any 9. In particular, 
Mhife) > 0. Let 



Then Mh{wg) = m for all 6. Moreover, it is easy to prove that Hh{w0) 
varies continuously from Hfi{g) to Hh{f) as 9 varies from to 1. This shows 
that for every E G (E'min(T?^, /i), -EmaxC^^, /i)), there is a function / with 
Mh{f) = m and Hh{f) = E. D 
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Lemma 9. Asn ^ oo, E^[j^{m, h, n) tends to E^i^{m, h) and E^ax{m, h, n) 
tends to Ema.xi^n,h). Moreover, in both cases, the convergence is uniform 
on compact subsets of [0, oo) (for the parameter m, keeping h fixed). The 
functions -Emaxl')^) and E^in{- , h) are absolutely continuous on [0, oo). 

Proof. Take any e > 0. Let / G C^ be a function such that Mfi{f) = m and 
Hh{f) < Emini'iTT'ih) + €. For each n, define a function /„ on Vn as simply 
the restriction of / on Vn. Then it is easy to see that M^^nifn) — ^ ^hif) 
and Hh^nifn) — ^ ^/i(/) as n — )• oo. Since e is arbitrary, this shows that 

limsup Emin{'nT',h,n) < £^min(m,/i). 

n— >oo 

Next, fix n and a function / E C^" such that Mh^nif) = "t- and Hh,n{f) < 
Emin{fn,h,n) + e. Define a function (7 on Z'^ as follows. For each x G V^, 
let fx be the translated function fx{v) '■= f{x + y), where the addition is 
modulo n in each coordinate. Let dVn denote the boundary of Vn when Vn is 
considered as a subset of Z*^. Recall the notation M^^nif^ U) from Section [5j 
Then 

^Mn,n{f.,dVn)=Yl J2hVix + y)\' 

xeVn xeVn y&dVn 

= Y, h''\fiz)\^\dVn\=Mh,n{f)\dVn\. 

Z&Vn 

This shows that there exists x £ Vn such that 

Mh,n{fx,dVn) < Mh,n{f)\dVn\n-'' < C{d,m)n-\ 
Take such an x and define 5 : Z*^ — t- C as 

if y G K, 




g{y) :- 

" otherwise. 

Clearly, M^ig) = m. Since x was chosen so that the 'boundary effect' 
is small, \Hh{g) — Hh,n{fx)\ < C{d,m,h)n~^. Since e is arbitrary and 
Hh,n{fx) = Hh,n{f) < Emin{m,h,n) + e, this shows that 

liminf £'min(r?2, /i, n) > E^[ii{m,h). 

n— ^-oo 

The proof for -Emax is similar. 

To prove uniform convergence on compact subsets of [0,oo), we will first 
prove it for compact subsets of (0,oo). We will show that the collection 
of functions {Emin{-,h,n))n>i is equi-Lipschitz continuous on any compact 
subinterval of (0, 00) and apply the Arzela-Ascoli theorem. This will also 
prove absolute continuity of the function -Emin on (0, 00). Continuity at zero 
follows from Lemma [8l The result for -Emax follows similarly. 

Take any n and < a < b. Take a < m! < m < b. Let / be an 
energy minimizing function (at grid size h) of mass m on the torus Vn- Let 
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/' := {m'/mf/'^f. Then Mh,n{f') = m' . Note that 

Thus, 

\HhAf) - HhAf)\ 

< \{m'/m) - l\GH,n{f) + lK/m)(P+i)/2 - l\Nh,M)- 
Since Nh^n{f) and Gh,n{f) can both be bounded above by 

this shows that 

-E'min("i-', h, n) < E^i^{m, h, n) + C{a, b,p, d, h)\m' — m\. 

Similarly, taking g' such that Mh^nig') = "m' and Hh,n{g') = Ejamiin' ,h,n), 
and letting 51 := {m/m')g', it follows that Eai\a{'ni-,h,n) is bounded above 
by ET^i^[m' ,h,n) + C{a,b,p,d,h)\m' —m\. This proves uniform convergence 
on compact subsets of (0,cx)). To prove uniform convergence on compact 
subsets of [0,00), one simply notices that the uniform bound on the energy 
given in Lemma [8] holds for functions on Vn as well (the bound will be same, 
independent of n) . D 

Recall that we say Uj ~ bj as j — )■ cxd if aj and bj are sequences that 
satisfy limj_j.oo o-j/^j = 1- 

Lemma 10. Let Cj be a sequence such that for some C > and a > 0, 
Cj ~ Cj'^ as j ^ CO. Then for any /3 > 0, 

fe-i „i 

V Cjik - jf ~ C7A;"+'^+^ / x°(l - xfdx as k ^ 00. 

Proof. Writing 

fc-i 



i=o 



fe-i 



j=0 



note that 



- E(c, A")(i - j/kf - - Y.{j/kr{i - j/k) 



j=0 

1^ Ic.--Cj"| 



j=o 



— I. 2-^ 



- I. Z^ 



^t^ a+1) 
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Since Cj/Cj" — )• 1 as j — )• oo, the above bound tends to zero as /c — )• oo. 
Riemann sum approximation gives 

lim ^ y (j7A;)°(1 - j/kf =C [ x"(l - xfdx. 

This completes the proof. D 

Lemma 11. Suppose that 1 < p < 1 + A/d. Then for any m > and A > 0, 

sup Emmim,h) < 0. 

0<h<A 

Proof. Fix m > 0. For each positive integer k, define a function /^ as follows. 
If |x|i > k, let fk{x) = 0. If |x|i < k, let 

fk{x) = Ak{k- |x|i), 

where A^ is a positive constant such that Mfi{fk) = m. Since the number of 
vertices x with |x|i = j asymptotes to C{d)j'^~^ as j — )■ oo, it follows from 
Lemma [TU] that 

Al ~ C(d, m)/i-'^A:-('^+2) as A; ^ oo. 

Note that if x and y are neighboring points, then fk{x) 7^ fk{y) if and only 
if one of them has i^ norm j and the other has l^ norm j ' + 1, for some 
< j < k. And in that case, 

\fk{x) - fk{y)\ = Ak. 

There are ~ C{d)j such pairs as j — )• 00. Thus, again by Lemma [T0| 

fc-i 
G{fk) ~ C{d)h^-'^AlY^j'^-^ ~ C{d,m)h'^-^Alk'^ - C((i,m)(/iA:)-2. 
i=o 
Again by Lemma [TOl 

fc-i 

^cid)h''Ai:^'^/-Hk-jr' 

j=0 
~ Cip, d, ^)/j<i-'^(P+l)/2fc-(<i+2)(p+l)/2^p+d+l 

Thus, for all k, 

Hifk) < Ci{p,d,m){hk)-^-C2{p,d,m){hk)-''^P-^^/\ 

If p < 1 + 4/d, then d{p — l)/2 < 2. It is now easy to see from the above 
bound that \{ h < A for some constant A, then 

^mm(m,/i) < iniH{fk) < -C{p,d,m,A) < 0, 

k 
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which concludes the proof. D 

Lemma 12. For any m > 0, £^jnax("i, h) = 2dm/ h? . 

Proof. Fix n. Let L, R and A be the matrices from Section [71 Let 

d 

Amax '■= max Ay = max > 4 sin (iryi/n). 
y&V„ yeVn -f— f 

Clearly, Amax < 4(i always, and Amax — )> 4d as n — >• oo. Note that for any 
/ E C^" with MhM) = ^^ 

ud-2 

HhM)<^^ E l/(^)-/(y)l' 

1^-2 

= ^{f,Tf) 

< ^-Amax 2^ |/(x)| = -^J^- 
x&V„ 

Thus, E^s,i^{m,h.,n) < 2dni/h?. 

Next, let / := vmhr^p^n/2]: where {py)y^v„ are the eigenfunctions defined 
in Section [7] and [n/2] is the element of Vn whose components are all equal 
to the integer part of n/2. Note that |/(x)| = vmh~'^n^''^ for each x. 
Therefore, Mfi^nif) = "t- ^'iid 

l,d-2 ud 

xr-jy 

= ^{f,Tf) - -^(m/i-^)(P+i)/2n-'^(P-i)/2 

Ld— 2 \ T d 

xeVn 

Since A[„/2] — ?" 4d as n — )• oo and the second term goes to zero, this shows 
that 

liminf E^maxCm-, /i,"-) > -nr- 

By Lemma m this completes the proof. D 



9. Large deviations for the gradient 

Recall the function ^^ defined in the statement of Theorem [5j The fol- 
lowing proposition summarizes some important properties of this function. 
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Proposition 13. The function ^d has the following properties: it is con- 
tinuous in {0,4d), it is strictly decreasing in {0,2d] and strictly increasing 
in [2d, id), ^d{2d) = 0, and 

lim^rf(a) = oo = lim ^^(a). 

In particular, ^^ is a continuous function from M into [0, oo] . 

Proof. It is obvious from the definition that ^ is strictly decreasing in (0, 2d) 
and by symmetry, strictly increasing in {2d, 4d). To extend the monotonicity 
up to the point 2d, we have to show that ^(a) > for all a < 2d. For each 
a <2d and < 7 < 1 define 



Kail) ■■-- 






so that 



(21) ^d(a) = sup Ka{j). 

0<7<1 

By dominated convergence, K^ is continuous on [0, 1) and differentiable in 
(0, 1), and has a right derivative at 0. Moreover, Ka{0) = 0. A simple 
computation shows that K'^{0) = 2d/a — 1 > if a < 2d. Thus if a < 2d, 
then Ka is strictly increasing at and therefore by (|2ip . ^^(a) is positive. 
By symmetry, ^^(a) is positive for a £ {2d, Ad) also. 

Next, let On be a sequence converging to a G {0,2d). By dominated 
convergence, Ka„ — )• Ka pointwise in [0, 1). It is easily seen that Ka„ is a 
strictly concave function for every n. By concavity and pointwise conver- 
gence, it follows that supg^^^^ -^0,1(7) converges to supo<^<i Ka{'y), proving 
that ^d is continuous at a. To show continuity at 2d, follow the same ar- 
gument and simply note that K2d attains its maximum at 0, since K2d is 
concave in [0, 1) and K2^{0) = 0. 

Lastly, note that for any fixed 7 E (0,1), Ka{'y) — )• 00 as a — )• 0. This 
proves that lima^^o ^d(a) = 00. By symmetry, limQ,f4rf ^d(a) = 00. D 

Next, let ^ = {C{x))x£V„ be a random function chosen uniformly from the 
unit sphere 



jneC^" : Y^ |'u(x)|2 = l|. 



xeVn 

The following theorem is a large deviation result for the gradient of ^ that 
is of fundamental importance in the rest of the manuscript. 



Theorem 14. For each a G {0,2d), 



lim ^logpf V |C(x) -C(y)|' < a ) = -*d 



(a). 
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Moreover, for any 6 > 0, the same limit holds for 



Vogpf Yl IC(x)-C(y)l'<«, max|e(x)|2<n-'^(i-^)Y 



^x,yeVn 



The same conclusions hold if a & {2d, oo) and the '< a' is replaced by '> a' 
in both expressions. 

Proof of Theorem^^ in the case < a < 2d. Fix a G (0, 2d) and a positive 
integer n. Let (p be the Gaussian random function defined in Section [H Let 
r, R and A be the matrices defined in Section [3 Let 

(22) r := R*^, 

where R* is the adjoint of R. Since i? is a unitary matrix, r has the same 
distribution as (p. Moreover, 

(23) Yl \<y)\' = E i'^(^)i' 

y&Vn xeV„ 

and 

(24) 5^ A,|r(y)p = (</>, r</<)= ^ \c^{x) - c^{y)\' . 

y£V„ x,yev„ 

x^y 

Let 

e(x) : = 



{T^yevMy)?) 



1/2- 



Then ^ is uniformly distributed on the unit sphere of CX" . Moreover, by 
(USD and dMD, 

^_ E-'2'6^" \(t>{x) -(/'(y)P 

E \ii^) - iiy)? = —^ — VTT2 

x,y&V„ 2^X(^Vn I'^^^'l 

_E,ey„A,|r(y)p 



Ej;ey„ l^(?/)l' 



Let ??(y) := n /i |r(?/)p. Then ry(y) is an exponential random variable with 
mean 1 and the ??(y)'s are independent. For any a G (0,2(i), 

( E i^(^) - ^(y)i' < «) = ip( E (^^ - ")i^(y)i' ^ « 

Y, (A, - a)ii{y) < o) . 






yeV„ 
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Thus, for any 9 G [0, 1/a), 



x,ySV, 



^l l + 9{Xy-a) 



yeVr, 
Note that we need the restriction that 6 < 1/a since Aq = 0. Now, 



yeVr 
n-1 

(25) 



n'^^'^U l + 0{\y-a) 



^ Y, \og(l-9a + WYs\n\7:y,/n)\=:In{9). 
yi,...,yri=0 ^ j=l 



The sequence of functions /« converges pointwise to the function / on 
[0, 1/a), where 

I (9) = [ log ( 1 - ^Q + 40 V sin^ (TTXi) ) dxi---dxd. 

This shows that for each 9 G [0, 1/a), 



(26) 



limsup^^iogpf V ie(^)-e(y)i'<«) <-^(^)- 






Now, 1(9) = Ka{9a), where Ka is defined in the proof of Proposition [T3l 
By dH]) and 1^, this shows that 

hmsupijlogpf ^ \^{x) - ^{y)\^ < a] 

(iyj\ ^x,yeVn, ' 

\^ i J x^y 

< - sup I (9) = - sup Ka{9a) = -^^(a). 

6»e(0,l/a) 0e(O,l/a!) 

This proves the upper bound in the case < a < 2d. Next, we establish the 
matching lower bound. 

As noted in Proposition [T3l Ka is a concave function and is strictly in- 
creasing at if a < 2d. Thus, the maximum of I in [0, 1/a) must be achieved 
either inside (0,1/a), or as —t- 1/a. The first case holds if I'(l/a) < 0, 
and the second happens if I'(l/a) > 0. The proof of the lower bound is 
different in the two cases. 

Case 1: /'(1/a) < 0. 

In this case, there is a unique 9* € (0, 1/a) where I is maximum. Fix 
any e > 0, and let 9' be a point so close to 9* from the right that I'{9') G 
(— e, 0). Fixing n, let r/' = ir]'iy))y£V„ be a cohection of independent random 
variables, where 77' (y) is an Exponential random variable with mean 1/(1 + 
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6'[\y — a)). Assume that rj' is defined on the same probabihty space as all 
other variables, and is independent of everything else. 

Recall the definition (f22]l of r. Let o"(y) := T(y)/|r(y)|. Since r(y) is 
a complex Gaussian random variable with mean zero, cr(y) and |T(y)| are 
independent random variables. Note that 



r{y) = a{y)\r{y)\=a{yy'''^y^ 



n' 



dh^- 



Define a random vector ^' = {C'{^))xeVn ^s follows: Take the variables a{y) 
and r]'{y) defined above, and let 



^■i^y-^i^^M 



Let (p' := Rt' , and let 



C'{x) :- 



1/2 /,r- |^/^„,^I2^V2■ 



{j:yevj4>'{y)\r'' {j:yev,y{y)\') 



Note that r', just like r, is a collection of independent complex Gaussian 

1 



random variables with E(T'(y)) = E(r'(y) ) = for each y, but unlike r, 



(28) E\T'{yt 



ndh'^{l + 9'{Xy-a))' 



Consequently, (p' is a complex Gaussian random vector satisfying 

but its components are not independent. 

Note that the map that takes {a, r]) to ^ is the same map as the one that 
takes {a, ij') to ^'. Therefore a simple change-of- measure computation shows 
that for any measurable set A C C^", 

P(eGA)=E(p(7?')lU'GA}), 

where 

Thus, if we fix 5 > and let E be the event 



E := |-4n'^e < ^ (A^ - aW{y) < 0, 

max|0'(x)p <n-°'(i-'^) ^ k'(2;)|4. 
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then 

Y. |e(^)-C(y)l'<«, max|e(x)|2<n-^(i-^)') 

^ ^ x,y&Vn ^ 

E[p(r?Ol|j;(A,-a)77'(y)<0, 



3^£ Vn I / 



(30) > np{v')lE) > e-'-'''^F{E) n YTm 



y^V.- ' -(^^-")' 



Now, if In is the function defined in ()25p . a simple computation gives 

1 ■^— \ Ay 

Thus, 



n'^J^^l + ^'(A,-a)- 



^{^,Y.^\-'^)^'(y))=in{e')- 



y&Vn 
But by independence, 

C{a,9') 



va^(;^E(^^-«)^'(y))^ 



n"' 



yeVn 

In particular, the variance tends to zero as n — )• oo. Since I'nid') — ^ -^'(^) ^ 
(— e, 0), this shows that as n — )• oo, 

(31) F(-4n'^€ < ^iXy- a)r]'{y) < J ^ 1. 

By the Gaussian nature of (/>' and (j29p . for any 5' > 0, 
lim P(nVmax|(/.'(x)p < n^'^) = 1, 

and similarly by (j28p and the independence of the coordinates of r', 

^ y n'^/i'^|r'(y)|2 ^ /" L^ dxi ■ ■ ■ dxd > 

in probability as n — )• oo. 
Combining the last two displays, and choosing 5' G (0,(5), we get 

(32) lim pfmax|0'(x)|2 < n-'^^^"^) V |r'(x)P^ = 1. 

n-s>oo \a;ey„ ^ / 



30 SOURAV CHATTERJEE 

Therefore, by ([30]), ([SlD and ([32]), F{E) ^ 1 as n ^ oo, and hence 

hminf^jlogPf y \C{x) - C{y)\^ < a, ma^\C(x)\'^ < n-'^^^-^n 

n-s>oo n" \ ^-^ x&Vn I 

^x,yeVn ^ 

> -Ae'e + hm i. V log(l + e'{\y - a)) 



y€Vn 



= -Ae'e-i{e'). 

Since e is arbitrary and 0' — )• 0* as e — )• 0, this shows that when a G (0, 2d) 
and /'(a) < 0, 

nminf47logP( y |C(x)-^(y)P <a, max |^(x)P < n^'^^^-^M 



xr^y 



> -I{0*) = -^d{a). 

Case 2: /'(1/a) > 0. 

Take any e > and let 6' solve 1 — 9' a = e. Fix n and define rj' , r', i^' and 
6 as in the previous case. Then as before, we arrive at the inequality (|30|) . 
with E being the same event. It suffices, as before, to show that 

(34) lim inf — r log ¥{E) > 0. 



n— >oo n 



Let 



and 



y(^Vn\{o} 



_ (Aq - a)V(O) _ aV(0) 

As before, it is easy to argue that Var(X„) — )• as n — )• oo. Again, 

A„ — a 



m.) = ^, E I 



'y 



. + d'(K-a) 
j/ey„\{0} ^ "^ y > 

— )• I' {9') as n — )• oo. 

Consequently, Xn — >• -?^'(^') in probability as n — )• oo. 

Next let ??"(0) be an independent copy of ?/'(0). Let r/" be the vector in 
C^" whose 0th component is r/"(0) and r]" {y) = r]'{y) for every y 7^ 0. Let 
r" and </)" be obtained from {a, r/") the same way as r' and (p' were obtained 
from {a,r]'). Then note that (p" is independent of 'r]'{0). Since the elements 
of the matrix R are bounded in absolute value by n~ ' ^ , and from definition 
we have 

(a) ^" = Rt", 0' = Rt', 

(b) r"(0) = a{0)^/r]"{0)/{nhY, t'(0) = a(0) Vr/'(0)/(n/i)^, and 

(c) r"(y) = r'iy) for all y / 0, 
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therefore for all x £ Vn, 

(35) |(/."(x) - ct>'{x)\ < n-'^/\nh)-'^/\'n"{0)y/^ - (7?'(0))i/2|. 

Fix 5' G (0, 6) and some e' so small that —d{l — 5') + e' < — d(l — 6), and 
define four events: 

Ei:={\Xr,-l'{9')\<e}, 

E2 := {-I'{d') - 3e < y„ < -I'{d') - e} 

< r? (0) < 

a a 

Es := /(n/i)'^max|(/)"(x)p < n^''^, 4 T" (nhyir" (x)\'^ > n'A, 

i^4:={|r?"(0)|<n'^}. 

Then by ([35]) . £^2 H £'3 n £'4 implies that for each x, 

\^'ix)\^<2\rix)\^ + 2\^'ix)-^"ix)\' 
<{nhy'^{2n^''^ + C{a,e)). 

Since t"{x) = t'{x) for all x / 0, this shows that for all n > C{a,e,6' ,e'), 
£^2 n E's n E4 implies 

max|,/.'(x)|2 < n-'^(i-^')+^' ^ |t'(x)|2 < n-'^^i-^) ^ \t'{x)\\ 

i'ev„\{o} xeVn 

Again, £^1 n £2 implies — 4e < X„ + 1^ < 0, which is the same as 

-4n'^e< ^(Aj,-aV(y)<0. 

yeVn 

Thus, for ah n > C{a,e,6' ,e' ,6), Ei D E2 D E3 D £4 implies E. So, to 
show (pl|) . it suffices to show that 



(36) Hm inf — r log p(£;i n £;2 n £;3 n £^4) > 0. 



n— >oo Ti"- 



Since r/'(0) is an Exponential random variable with mean 1/e, it is easy to 
see that 

2n<'e^ ( en''{I'{e') + 3e) 
™(£'2) > exp 



a \ a 



As argued above, X„ — )• I'{6') in probability and therefore P(£^i) — )■ 1 as 
n — 7- 00. The probability of £4 tends to 1 trivially, and P(£?3) — )■ 1 by the 
same logic that led to ([32]). Thus, F{Ei n £^3 n £'4) -^ 1. Lastly, observe 
that the events Ei, E^ and S4 are jointly independent of £^2- Combining all 
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of these observations, we get 



lim inf — r log F{Ei n £^2 n S3 n -£4) 

n— >oo n 

= liminf -T(logP(Si n £^3 n £4) + logP(£;2)) 

n— >oo n 

e{r{e') + 3e) 
a 



> 



Since e is arbitrary and < I'{0) < I'{0) = 2d— a for all 9 (by concavity and 
the assumption that /'(1/a) > 0), this proves ([36]) and hence ([M|) . leading 
to the proof of p3]) when /'(1/a) > 0. Combining ([27|) and (f33]) . the proof 
of Theorem 1141 for a E (0, 2d) is complete. D 

Proof of Theorem^^ in the case 2d < a < 00. The proof is very similar to 
the previous case, with a few important modifications. Let all notation be 
as before. Note that for any a G {2d, Ad), 



,ysv„ ^ ^yGV„ 

Thus, for any (9 G [0, l/(4(i - a)). 






J-J: 1-^fA, 



The condition 6 < l/(4(i — a) ensures that the right-hand side makes sense, 
since Ay is uniformly bounded above by 4d. Now, 



^'^°'n l-^(A,-a) 



n— 1 / d 



-^ E logfl + 0a-40E^^"'(^2/i/n)j =: Jn(e), 

«i,...,«rf=0 ^ i=l ' 



and the sequence of functions J„ converges pointwise to the function J on 

[0, l/(4(i- a)), where 

J{0) = log ( 1 + 6^0 - 46* E sin^ (^2;j) ) dxi • • • dx^. 



/ logf 1 + 6*0 -46* Vsin2(7rxi)) 
■^[0,1]'* V ^ / 
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This shows that for each 6 G [0, l/(4(i — q)), 

hmsup^logpf Y, I^W-^(y)|'>a) <-J{9), 

and therefore 

Hmsup^logpf V Ux) - i{y)\^ > a] < - sup J{e). 

n^oo n K^^y^vn, ^ 9e(0,l/(4d-a)) 

Now, putting 7 = (Ad — a)9 gives 

sup J{9) = sup J(7/(4(i — a)) 

6>e(0,l/(4d-a)) 7e(0,l) 



e(o,i)J[o,i]^ V 4d-a4d-a^ y 



= sup 

7e{o,i) 

The above expression is the same as ^rf(4(i — a), except that we have cos 
instead of sin. However, this does not matter, since we can break up the 
hypercube [0, 1]*^ into a union of smaher hypercubes hke [ai, ai + 1/2] x • • • x 
[ad,ad + 1/2] where each Oj G {0,1/2}, and then, within each hypercube, 
replace cos by sin by a change of variable yi = 1/2 — Xi when Cj = and 
yi = 3/2 — Xi when Oj = 1/2. Thus, 



(37) 



limsup^logpf V \^{x) - ^{y)f > a) 

n—>-oo n \ '--' J 






< - sup J{Q) = -^di^d - a) = -"^dioi). 

6»G(0,l/(4rf-a)) 



Next, we turn our attention to the lower tail. The function J is continuous 
in the interval [0, l/(4d — a)) and differentiable in (0, \/{Ad — a)). Moreover, 
J is differentiable from the right at and differentiable from the left at 
l/(4(i — a), and an easy computation gives J'(0) = a — 2d. Since a > 2d 
and J is continuous at 0, this shows that J must be strictly increasing 
in a neighborhood of 0. It is easy to check that J is a concave function. 
Since J is increasing at 0, its maximum in [0, l/(4(i — a)) must be achieved 
either inside (0, l/(4(i — a)), or as ^ — )■ l/(4(i — a). The first case holds if 
J'{\/{Ad — a)) < 0, and the second happens if J'(l/(4(i — a)) > 0. Just as 
before, the proof of the lower bound is different in the two cases. 

Case 1: J'(l/(4d-a)) < 0. 

In this case, there is a unique 6* G (0, l/(4(i — a)) where J is maximum. 
Fix any e > 0, and let 0' be a point so close to 9* from the right that 
J'i9') G (—6,0). Given n, let (?7'(2/))y6V„ be a collection of independent 
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random variables, where ri'{y) has the Exponential distribution with mean 
1/(1 — 9'{Xy — a)). (Note the minus sign in front of 6', which was plus in the 
case a < 2d.) As before, assume that rj' is defined on the same probability 
space as all other variables, and is independent of everything else. 

Given r/', define r', (p' and ^' as before. Then all the properties of these 
vectors are same as before, except that now 



E\T'{y)\ 



{nh)d{l - e'{\y - a)) 



and 



Defining 



^\^'^^^\'^ inkYii-e'iM-.)) fo-^^--^^»- 



J-i l-e'{\y-a) 



Fix (5 > and let E be the event 



Then as before we arrive at the inequality 

E le(x)-e(y)P>a, max|e(x)|2<n-'^(i-'^)) 



max|,/.'(x)|2<n-'^(i-^) 






(38) > e-4^'"''T(i?) n Y^\ 



y^V.' -(^^-«) 



Exactly as before, we can now argue that 

—1 2, {\ ~ '^)v'{y) ~^ ~J'{Q') ill probability as n — )■ oo. 






Since J' {6) G (— e, 0), this shows that as n — t- oo, 

(39) P( < ^{\y- a)r]'{y) < 4n"'e J -^ 1. 
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The inequality (f32]l continues to be valid, and therefore, by (f38]) . (|39]) 
and (1321). 



lim^inf ^logpf Y^ \Cix) - C{y)\'^ > a, max\^{x)\^ < n''^^^-^n 
> -4e + lim ^ V logfl - e'(Xy - a)) 

n— >oo n ^ — ' 



yeVn 



= -4e-J{e'). 

Since e is arbitrary and 0' — t- 0* as e — t- 0, this shows (as in (j37p ) that when 
a £ (2d,4d) and J' {1/(4.(1 - a)) < 0, 

liminfi^logPf y \^{x) - ^{y)\^ > a, max |e(x)|2 < n-'^^^-^M 

(40) "^~ ^ v.fit„ -ey„ y 

Case 2: J'(l/(4(i - a)) > 0. 

Take any e > and let 6' solve 1 — 0'(4(i — a) = e. Fix n and define 
rj' , t' and ^' as in the previous case. Let [n/2] is the vector in C " whose 
components are all equal to the integer part of n/2. Then as before, we 
arrive at the inequality psp . with E being the same event. Let 



and 



^^'■=7;! Yl ^^y - «)^'(y)' 



yeVn\{[n/2]} 



(A[„/2]-a)7?'([n/2]) 



Y 

n"- 



Note that \nj2] — )• 4fi as n — )• oo. 

As before, it is easy to argue that \ai{Xn) — ?■ as n — )• oo. Again, 



nd Z^ 1 _ 0'fx - a) 
-J'(9') 



as n — )• oo. 



Consequently, Xn — )• —J'{0') in probability as n — ;■ oo. 

Next let T/"([n/2]) be an independent copy of 77' ([n/2]). Let r/" be the 
vector in C " whose [n/2]th component is 7]"{[n/2]) and r]"{y) = r]'{y) for 
every y ^ [n/2]. Let t" and cp" be obtained from (cr, 77") the same way as 
r' and (f)' were obtained from (cr, 77'). Then note that (p" is independent of 
r/'(0). Exactly as we proved (f35]l . it follows that for all x £ Vn, 

10" (x) - 0'(x)| < 7i-'^/2(n/i)-'^/2|(7?"([n/2]))i/2 - (?7'([n/2]))i/2|. 
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Fix 5' G (0, 5) and some e' so small that —d{l — 5') + e' < —d{l — (5), and 
define four events: 

Ei:={\Xn + J'{e')\<e}, 

E2 := {J'ie') +€<Yn< J'{e') + 3e} 

\n/2] - " Hn/2] " " 

^3 := |(n/i)Vax|0"(x)|2 < n^''^, \ ^ (n/i)'^|r"(x)p > n"^'), 

I X t V Ti ''' I 

•^ xeV„\{[«/2]} ^ 

E4:=K([n/2])<n'^}. 

Then £'2 H E's Pi -E4 implies that for each x, 

< (n/i)'^(2n^''^ + C(a,e)). 

Since t"{x) = t'{x) for all x 7^ [^/2], this shows that for n > C{a,e,5' ,e'), 
E'2 n E'3 n E4 implies 



max \<P'{x)\^ < n-^(i-^) Yl l^'(^)l' ^ ^"''^'"'^^ J^ |t'(x) 



a;6V„\{0} aieVn 

Again, £^1 n -E'2 implies < Xn + ^ < 4e, which is the same as 

< J] (Aj, - aWiy) < Ve. 

J/6V„ 

Thus, for n > C{a, e, 5' , e', (5), £'1 n E'2 n £"3 n E4 implies E. So it suffices to 
find a lower bound for the probability of Ei D E2 D E^ D E4. 
Since r]'{[n/2\) is an Exponential random variable with mean 

1 



l-6''(A[„/2] -a)' 

and A[„/2] — > 4(i as n — )■ 00 and by definition of 6', 1 — 0'{4:d — a) = e, it is 
easy to see that 

,i„,int^l„gP(E,)>-im±M^ 

n-5>oo n"- Ad — a 

The proof is now completed exactly as for the lower tail in the case a G 
(0,2(i). D 

The following theorem is the main result of this section. 

Theorem 15. For each e > and a > 0, 

|2 



lim —7 log 

n— ^00 fi'* 



E i^(^)-^(2/)i' 



a 






< e = -*d,,(a). 
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where 



*d,.(a) 



^rf(a + e) ifa<2d-e, 
^rf(a - e) ifa>2d + e, 
if2d-e<a<2d + e. 



In particular, "$^,6 converges uniformly to ^^ on compact subsets of (0, Ad) as 
e — 7- 0. As in Theorem WAi the same limit holds if we include the additional 
requirement that vn&yix^Vn IC(^)P — ri~'^^^~ ' . 

Proof. Obvious from Theorem [T3] and Proposition [T3l D 

10. The variational problem 
For each £' > 0, ?n, > and /i > define 

'2h'^E 



Q{E, m, h) := log m — ^a 



m 



When 7TT- = 0, the right-hand side is interpreted as — oo. With this definition, 
it is easy to verify that for fixed /i, is a continuous function from [0, oo) x 
[0,oo) into [—00,00). 

Given tuq > and i?inm("i, /i) < Eq < £'max('7i, /i) let M.{Eo,mo,h) 
denote the set of all {E,m) that maximize Q{E,m,h) in the set 

TZ{EQ,mo, h) = {{E,m) : < m < mo, 

max{E~{m,h),0} < E < E+{m,h)}, 



(41) 
where 

(42) 

Define 



E (m,/i) := E'o --Emaxlm-o -m, /i), 
E^{m, h) := Eq - Ejamimo -m,h). 



@{EQ,mo,h) := max @{E,m,h). 

(E,m)e'R,{Eo,mo,h) 



The following lemma lists some important properties of the sets TZ and A4. 

Lemma 16. Suppose that -Emin("i-0) ^) 1^ Eq < dmo/h'^. Then the set 
TZ{EQ,mQ,h) is a non-empty compact subset of IB?, and so is Ai{EQ,mQ,h). 
Moreover, any (E,m) £ A4{EQ,mQ,h) satisfies m G (0, mo] CLnd E = 
E^{m,h), where E~^ is defined in (j42p above. 

Proof. From Lemma[8]and Lemma[9]it follows that TZ{Eq, mo, h) is simply a 
region enclosed between the graphs of two continuous functions on a closed 
interval and therefore, is a compact subset of M^. It is clearly non-empty. 
By continuity of 0, this shows that M{EQ,mQ,h) is also compact and non- 
empty. 

It is obvious that m > 0, since Q{E,m,h) = —00 if m = 0. Let E~ and 
E'^ be as in (I42p . Let E*{m,h) := dm/h"^. Lemma [T2] implies that for any 
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m £ [0,mo], 

drriQ 2d{mQ — m) 



< 



/l2 /i2 

dm d{mQ — m) 



< E*{m,h). 



/i2 /i2 

Moreover, E* (m, h) is non-negative, and hence 

E*{m,h) > max{E~{m,h),0}. 

Now fix any m G [0, mo]. If E^{m, h) < 0, tlien tfiere exists no E sucli that 
(E, m) G TZiEo, tuq, h), and hence no E such that (E, m) G TZ{Eq, tuq, h). So 
assume that E~^{m,h) > 0. By Proposition I13[ Q{E,m,h) increases strictly 
as E increases from to E*{m, h), and then starts decreasing as E increases 
further. Therefore, if we impose the restriction that [E,m) G TZ{EQ,mQ,h), 
then for fixed m, Q{E,m,h) is maximized at in.m{E^{m,h),E*{m,h)}. 

Suppose that (£', m) is a point in A^(ii^O) w-O; h) such that E = E*{m, h) < 
E~^{m, h). This is clearly not true if m = mQ, since E*{mQ, h) = dniQ/h? > 
Eq = E^{mQ, h). We claim that this is impossible even if m < rriQ. Indeed, 
if this is true for some m < ?n,o, then since -Emin is a continuous function 
by Lemma m we can choose a slightly larger m! > m such that E*{m' , h) < 
E+{m',h). But then E*{m',h) G [E-{m' ,h), E+{m' ,h)], and 

e{E*{m',h),m',h)=logm' -^i'J^^^-^^^^T^ 

\ m 

= \ogm' > logm, = Q{E,m, h), 
showing that {E,m) cannot belong to M{Eo,mo,h). D 

11. Upper bound 

Fix /i > 0, and some Sq £ 1^ and mo > 0. The numbers p, d, h, Eq 
and mo will be fixed throughout this section and will be called the 'fixed 
parameters'. Any constant that depends only on the fixed parameters will be 
denoted simply by C, instead of C{p, d, h, EQ,m.o). If the constant depends 
on additional parameters a,b, . . ., then it will be denoted by C{a, b, . . .). 

Recall the random function cp defined in Section [6] and the objects M., TZ, 
@ and Q defined in Section [TOl 

Take any positive integer n. For any 6 > 0, we will call a function / G C^" 
a (5-soliton if there exists a g £ C^" such that 

(a) IIZ-s-lloo < S, and 

(b) there exists {E*,m*) G M.{EQ,m,Q,h) such that 

\{Eo - E*) - Hh,n{9)\ < S and 
\{mo - m*) - Mh^r^{g)\ <S. 
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Theorem 17. For arbitrary e, 5 G (0, 1), let B = B{€,6,n) be the event 

{\Hh,ni4>) ~ -E'ol < e, \Mh^ni(p) — "T-ol ^ £> 4> '^^ ^ot a (5-soliton}. 

Then for any fixed 6 £ (0, 1), 

logF(B(€,6,n)) ^^^ 

lim sup lim sup ^ < 1 — mo + fc)(i?o, w-Oi ")• 

e— >-0 n—^oo ri 

The strict inequahty is the main point of the above theorem. Let us now 
embark on the proof of Theorem [T71 We begin with two simple technical 
lemmas. 

Lemma 18. Let zi, . . . , z^ be standard complex Gaussian random, variables, 
where k > 2. Let S := X]j=i k«P- Then for any x > 2 and < y < x/2, 

(43) P(|S'-x| <y) = e:^p{k + klog{x/k) - x + R{x,y,k)), 

where 

Ckv 
\Rix,y,k)\ < Clogk + Clogx + - + y, 

X 

where C is a universal constant. Moreover, 

(44) V{S <x)< exp{k + k log{x/k)). 

Proof. The random variable S has a Gamma density with parameters k 
and 1. Explicitly, the density function of 5 is: 

P^'^ = Jk^Ty.- 

Let T(x, k) := k + k \og{x/k) — x. Note that for any t G [x — y,x + y], 
\\og p{t)-T{x,k)\ 

< I log{k — 1)! — (/clog A; — A;)| + | logt| + k\ logt — logx| + |x — t| 

Ckv 

< C log /c + C log X H h y. 

X 

(The assumption that y < x/2 was used to bound second and third terms.) 
Further, note that 

fx+y 



rx+y 

P(|5-x| <y)= I p{t)dt 

J x—v 



ix-y 

^^T(x,k) r '\log p(t)-T(x,k)^^_ 
Jx-y 

Using the bound from the previous display finishes the proof of the first part 
of the lemma. To prove the second part, note that 

fX .k-l -t 

- Jo {k-l)l^^ " ¥ - k^^' 
This completes the proof. D 



40 SOURAV CHATTERJEE 

Lemma 19. Let "$([,€ be the function defined in Theorem\T5i Then for any 
a G [0, oo) and any L > 0, 

*d,e(a) > min{*rf(a), L} - a(e, L, d), 

where a(e, L, d) is a quantity that depends only on e, L and d (and not on 
a), such that for any fixed L > 0, lim^^^o o(£) L, d) = 0. 

Proof. In this proof, a(e, L, d) will denote any constant with the properties 
described above. 

By the properties of ^^ listed in Proposition [131 there exist < ci < 
C2 < 4d such that ^d(a) > L whenever a [ci, C2]. Fix c'l and c'2 such that 
< c'l < ci < C2 < c'2 < 4:d. Note that ci, C2, c'l and C2 can be chosen 
depending only on L and d. By the uniform convergence of ^d,e on compact 
sets, if a G [c'l, c'2], 

^d,e(a) > ^d(a)-a(e,L,d). 

On the other hand, by the definition of ^d,€j it is easy to see that if a 
[c[, C2], for all sufficiently small e (depending only on L and d), '^d,eict) > L. 
In particular, if a [ci,C2], then 

^rf,,(a) >L-a{e,L,d). 

The proof is completed by combining the two cases. D 

For a subset U of Vn, recall that If^ denotes the set Vn\U, and dU denotes 
the set of all vertices in U that are adjacent to some vertex in U^. Recall 
also the definitions of Mh,n{f,U), H^M^U), N^M^U) and GuM^U) 
from Section [8l The following lemma shows that the nonlinear component 
of any / must come from a small region. 

Lemma 20. Take any f G C ". Let m := M^^nif)- Then, given any 
€ G (0, 1), there exists a non-empty subset UofVn such that 

(i) 1^1 < ]^- 

(ii) |/(2;)| <efor all x e U". 

(iii) Nh,n{f,U^)<eP-'m. 

(iv) Mh,n{f,dUUdU'') < 2em. 

Proof. Let Ui be the subset of V on which / is bigger than e. If this set is 
empty, let Ui be any singleton subset of Vn- Then note that / < e outside 
Ui, and 



Also note that 



l^^ll — TTo - 



Define U2,U^,... as follows: For each i >2, let Ui be the set of vertices that 
are either in Ui-i or adjacent to a vertex in C/j_i. Note that C/j_i C [/j for 
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each i. Let Wi := dUi. Then Wi,W2, ■ ■ ■ are disjoint sets, and Wj+i = dU^ 
for each i > 1. Thus, for any k, 

l<t<k ^ — ' l<i<fc ■^ — ' 

xedUiUdU^ x€WiUWi+i 

2h'^ v-^ v-^ , ^ , ,,9 2m 

^— EE 1/(^)1^ X- 

Since each element of Uk is within i^ distance k — 1 from some element of 
Ui, and the i^ ball of radius k — 1 around any vertex has < {2k — 1) points, 

max \Ui\ = \Uk\ < (2k - l)'^\Ui\. 
i<i<k 

Choose an integer k such that 1/e < k < 2/e. The proof is completed 
by choosing an i between 1 and k that minimizes h'^J^x&dUudU" I/^^I^j ^'^d 
defining U to be the set Ui. D 

Given any / S C^", recall that U{f,e) denote the set of all non-empty 
subsets of Vn that satisfy conditions (i) through (iv) of Lemma [20l Then 
Lemma [20] savs that V({f, e) is non-empty for any e > 0. 

Lemma 21. For any e G (0,1), E >0 and m, > 0, let K = K{e,n,E,m) 
be the event 



{\Mh^n{(p) — "T-ol ^ £) \Hh,ni4') ~ -^ol ^ £) aiid for some U &U{ 
|M^,„(</), C/^) - m| < e, |G/,,„(0, f/^) - E\ < e}. 
T/ien /or any L > 0, 

logP(i^(e,n,£^,m)) 



^5 "-y ; 



lim sup ■ 



n^ 



< max{l — mo + 6(-E, m,, /i), C — L} + a(e, L), 

where C depends only on the fixed parameters, and a(e, L) is a quantity that 
depends only on e, L and the fixed parameters, such that for any L > 0, 
lime_!.o a(e, -L) = 0. In particular, a{e,L) does not depend on E and m. 

Proof. Fix L>0, -E>0, m>0 and e G (0,1). Throughout this proof, 
a(e, L) will denote any constant with the properties outlined in the statement 
of the lemma, and o(l) will denote any constant that depends only on e, E, 
m, L, n and the fixed parameters, that goes to zero as n — )• oo while keeping 
the other parameters fixed. 

Choose a positive integer n and a set U '^Vn- For notational simplicity, 
define 

M := Mh,n{^,U), 

Hh,ni^,U), 
Mh,n{(l>,Ul, 



G' 
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Let (j)' be an independent copy of </>, and define 

, , \(t)'{x) if X G [/, 
t[x) := < , , 
^ ' \^{x) ifxGC/^ 

Note that r has the same distribution as (f), and is independent of {(l){x))x<^u- 
Let Kq = KQ{e, n, U, E, m) be the event 

{\Mh,n{(l)) - mol < e, \Hh,n{<P) - ^o| < e, |M' - m| < e, 

|G'-^| <e, t/G^/((/),e)}. 

Note that it is possible for Kq to happen only if m < M' + e < Mh^ni^P) + e ^ 
mo + 2e < C and E < G' + e < CM' + e < C. Therefore we will assume 
these upper bounds on E and m in what follows. We will also assume that 

(46) \U\ < -p^dT^' 

since without this condition, the event Kq is impossible. Let Ki be the event 

2\U\ 



Ki := Mh^r^iT, U) < 



.d 



Since M + M' = Mh^n{4'), therefore if Kq happens, then 

(47) \M - {mo - m)\ < \M + M' - mo| + \M' - m\ < 2e. 

Again, if Kq n Ki happens, then 

\Mh,n{T) -m\< \M' -m\+ Mh,n{r, U) 
2\U\ 



< e + 



n^ 



(48) <e + Ce~^'^+^^n~'^=:ei, 

and 

\Gh,n{^) -E\< \G' -E\ + G;,,„(r, U) + ^ ^ 



2 



W -r P 



xeU xedU" 

(49) <Ce + C7e-('^+2)n-'^=:e2. 
Let 

pi := P(|M - (mo - m)\ < 2e), 

P2 := n\Mh,n{T) -m\< ei, \Gh,n{r) - E\ < €2). 

Note that Kq and i^i are independent events, and there is a positive uni- 
versal constant Cq such that F{Ki) > 1/Cq. Thus, by (gT]), (jMD, (USD and 
the independence of M and r, 

(50) P(Ko) < CoF{Kq)F{Ki) = CqF{Kq n Ki) < CqpiP2. 



INVARIANT MEASURES AND THE SOLITON RESOLUTION CONJECTURE 43 

Define 



e(x) := 



TiX 



{EyevMy)\'f'' 
Then ^ is uniformly distributed on the unit sphere of C^" . Note that 



2h'-''Gh,n{i) = E l^(^)-^(y)l' 



xr^y 



Mh,n{r) 
2h^Gh,n{r) 



E \T{x)-T{y)\' 



x,y£Vn 



Mh,n{T) 

Thus, if |M/i^„(r) — m| < ei and |G/i^„(r) — E\ < €2, then (since m < C and 
S < C, as observed before), 



2h'-''Gh,n{0 



2h^E 



m 



2h' 



Gh,n{^) 



E 



Mh,n{T) m 

^ Cei + Ce2 

< -■ ^ =: £3, 



{m — ei)m 

provided that m > ei. If ?tt, < ei, define €3 = 00. Now, it is a simple 
probabilistic fact that ^ and Mh^ni^) are independent. Hence p2 < PzPa-, 
where 

P3:=P(|iV4,„(T)-m|<ei), 

2h'^E 

Thus, from (jSOh we have 

(51) 



2h'~^GhAi) 



m 



<e3 • 



P(i^o) < CoPlP3P4. 



Our next task is to get upper bounds for pi, ^3 and ^4. To bound pi, we 
consider two cases. First, if mo— w- > 2e, then we apply (j43p from Lemma [T8l 
(with k < C(e), x = n'^{mQ — m) and y = n'^e) to get 



logpi 



n" 



-(mo - m) + e + o(l) 



where recall that the notation o(l) stands for a quantity depending only on 
e, E, m, L, n and the fixed parameters, that goes to zero as n — )• 00 with 
all else fixed. In particular, o(l) does not depend on our choice of U. 

Next, if mo — rn < 2€, apply ()44p from Lemma [18] (with k < C(e) and 
X = Sn'^e) to get 

logpi 



ol . 



n" 



Combining the two cases, we get 
logpi 



(52) 



n" 



< — (m-o — m) + 3e + o(l). 
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We deal with ps similarly, lim > 2e^'^, we apply (H3]l with k = n'^, x = n'^m 
and y = n e\ to get 

(53) i^ = l + logm-m + Ce3/4 + o(l). 
When m < 2e^'^, we apply ()44p with k = n and x = 6n e to get 

(54) l^<i + C7ioge + o(l). 

Again, if ?Ti > 2e^''*, then £3 < Ce^" + o(l). Therefore by Theorem 1151 we 
have that if m > 2e^'^, then 

logp4 . ,., /2/i2E\ 

By Lemma [121 this gives 

(55) i^<max|-M/,(^), -l} + a(6, L) + o(l), 

where recall that a(e, L) stands for a quantity that depends only on e, L 
and the fixed parameters, that goes to zero as e — t- for any fixed L. In 
particular, a(e, L) does not depend on E, m, n or U. 

Combining ([5T]) . ([52]) . ([53]) and ([55]) and the observation that m < C, we 
see that when m > 2e^'^, we have 

logF(Ko) ^ logpi + logp3 + logp4 
n*^ — ^d 

< 1 — niQ + log m 

+ maxi -^rf ( j , -l\ + a(e, L) + o(l) 

< max{l - mo + e(^, m, /i), C - L} + a(e, L) + o(l). 
On the other hand, since 

Cloge < -L + a(e,L), 

and logp4 < and m < C, it follows from ([5T]) . ([52]) and ([M]) that when 
m < 2ei/'^, 

logF(Ko) ^ logpi + logp3 

^d — ^d 

< C-L + a(e,L) + o(l). 

Combining the last two displays, we see that for all e G (0, 1), ra > 1, i? > 0, 
m > and [/ satisfying ([16]) . we have 

logF(Ko(e,n,[/,g,m)) 

< max{l - mo + 6(£;, m, /i), C - L} + a(e, L) + o(l). 
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Now note that K can be written simply as 

K= \J Ko{e,n,U,E,m). 

ucv„ 

Since there are at most e ''^•' ^^^ terms in the above union, this completes 
the proof of the lemma. D 

Lemma 22. Fix n and let K = K{e,n,E,m) be the event defined in 
Lemma [21]. // K happens, then there exists a function rj on Vn such that 

(a) \\4>-v\\oo < e, 

(b) \{Eo - E) - Hh,n{v)\ <2e, and 

(c) \{mo-m)-Mh,n{v)\<Ce + C€P-\ 

Proof. Suppose that K has happened. Choose any U G U{4>,e) satisfying 
the conditions of K and let M, H, M' and G' be as in (j45p . Let rj be the 
(random) function 

if X G [/, 
if a; G V. 



max \(b{x)\ < e. 




M — (mo — m)\ 

\Mh,n{(t)) - M' - {mo - m)\ < 2e, 



\Hh,n{v) - (Eo - E)\ < \Hh,n{<P) - Eo\ + \G' - E\ 

(58) <2e + C Y, l^x)]"^ + C Y, \^(^)\^'^^ 

<C7e + CeP-\ 

This completes the proof. D 

Lemma 23. Let K{e,n,E,m) be as in Lemma ^2T[ Given any mi > and 
El > 0, there exists Go = Go{Ei,mi) and Gi such that if m > mi, E > Ei, 
e < Co and {E,m) is at i°° distance greater than Gie + Cie^"^ from the 
set TZ{Eo,mo,h), then for all n > G2{Ei,mi,e), the event K{e,n,E,m) is 
impossible. 
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Proof. If K{e,n,E,m) happens, then by Lemma [22| there exists a function 
r] eC^" such that 

(59) \{mo - m) - Mh,^{r])\ < 2e, 

(60) I (So -E)- Hh,n{v)\ <Ce + Ce^'K 

Let E' := EQ — Hh^nif]) s-iid m! := mo — Mh,n{v)- Then the above inequahties 
may be rewritten as 

(61) \m-m\<2e, \E - E'\ < Ce + Ce^-^. 

Now, if ?n, > mi and if e is small enough depending only on mi and the 
fixed parameters, then by (f59l) . 

(62) m' = mo - M,,,„(r/) > m - 2e > mi - 2e > 0. 

But by definition, m' < mo- Therefore, m' G [0, m,o]. Next, note that 

Eram{Mh^ri{rO,h,n) < Hh,n{v) < -E^max(M;,_„(r?), /l, n), 

which is the same as 

Eq - £'max("T'0 " "^', /i, n) < £^' < E'o - £'min(mo - m', /i, n). 

Again, if £^ > Ei, then for sufficiently small e (depending only on m,i and 
the fixed parameters), (IGOp implies that 

(63) E' = Eo- Hh,n{v) >Ei-Ce- Ce^'^ > 0. 

Combining the last two displays, we see that ii E > Ei, m, > mi and e is 
sufficiently small, then 

maxjE'o - -E'max(m'0 - m', /i, n), 0} 
(64) 

< E' < Eq- £'min("10 - m', /l, Tl). 

By Lemma [9] and the characterization (|4ip of the set TZ{Eo,mQ,h) from 
Section [lOl the equations (f62]) . ([63|) and (fM|) show that if e < Co(-Ei,m,i) 
and n > C2(£'i,m.i, e), then the event K{e,n,E,m) implies that the point 
{E' ,m!) is within distance e from the set TZ{EQ,m.Q,h). By (f6T]l . the proof 
is done. D 

Lemma 24. For any e G (0,1) and any closed set A C [0, oo)^, let F = 
F{e,n,A) be the event 

{\Mh-n{(p) - mo\ < e, |-fffe,n(0) - -E'ol < e, and for some U £ U{(j), e). 

Then for any such set A, 

limsuplimsup j — ^ < 1 — m-o + @{EQ,mo, h,A), 

e— >0 n— >oo n, 

where 

Q{Eo,m,o,h,A) := sup Q{E,m,h). 

(E,m)eAnn(Eo,mo,h) 
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Proof. If F{e, n, A) happens, then M;j^„((/), V^) and Gh,ri[4>^ U*^) are bounded 
by constants depending only on the fixed parameters. Also, TZ{Eq, mo, h) is 
a bounded set. Hence we can assume without loss of generality that A is 
contained in some bounded region determined by the fixed parameters. 

Then A is compact, and therefore there is a minimal collection C(e) of 
points in A such that the union of ^°°-balls of radius e around these points 
covers ^4, and the size of this collection is bounded by a constant depending 
only on e and the fixed parameters. Then we have 

F(e,n,A)C (J K{e,n,E,m), 

(E,m)£C{e) 

and therefore 

logP(F(e,n,A)) log|C(e)| ^ log¥ {K{e,n,E,m)) 

(00) -j j^ "5 ~r II13,X -J . 

W^ n'^ (£;,m)ec(e) n*^ 

Fix L > 0. Fix £^i > and m\ > so small that whenever E < E\ ox 
m < mi, 

(66) l-mo + Q{E,m,h) <-L. 

(It is easy to see that this is possible, by first choosing m-i so small that 
1 — rriQ + logm-i < — L, and then choosing Ei depending on mi.) Let 
Co = Co{Ei,m.i) be the constant from Lemma [23l and assume that e < Cq. 
Let Ci be the second constant from Lemma [23j Let C'(e) be the set of 
points {E,m) G C(e) that are at i°° distance < Cie + Cie^"^ from the set 
Tl{EQ,m.Q,h), and satisfy E > Ei and m > m.i. Let C"(e) be the set of 
points in {E,m) G C(e) such that E < Ei oi m < mi. Let C"'{e) be the set 
of all remaining points in C(e). Then by Lemma[23l for each {E, m,) G C'"(e), 

logF( K(e,n,E,m)) 
lim -^ = — oo. 

For any {E,m) G C"(e), Lemma[2T]and the inequality (j66p imply that 
logP(i^(e,n,S,m)) 



lim sup ■ 



j^d 



< max{l — TTio + max Q{E,m,,h), C — L} + a{e,L) 

{E,m)eC"{e) 

<C-L + a{e,L), 

where a(e, L) is a quantity depending only on e, L and the fixed parameters, 
such that for each fixed L > 0, lime_j.o a(e, L) = 0. 
Finally, again by Lemma [211 for any {E,m,) G C'(e), 

logP(K(e,n,£;,m)) 
lim sup ^ 

n— >oo ?^ 

< max{l — TTT-o + max Q(E,m,h), C — L} -\- a(e,L). 

{E,m)€C'{e) 
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Combining the last three displays with (I65p . we get 

logP(F(e,n,A)) 
hm sup 2 

n—^oo ^ 

< niaxjl — mo + max Q(E,m,h), C — L} -\- a(e,L). 

{E,m)eC'{€) 

Now, by the definition of C'(e) and the continuity of B, it follows easily that 



Therefore, 



limsup max Q(E,m,h) < Q(E,m,h,A). 

,^0 {S,m)GC'{.) 



logP(F(e,n,^)) 
lim sup iim sup -j 



>0 n—>-oo IT' 

< max{l - mo + Q{E,m,h, A), C - L}. 
Since L was arbitrary, this completes the proof. D 

Lemma 25. Take any e, 5 € (0, 1). Let As be the set 

{iE,m) £ [0,oo)2 : max{\E - E*\,\m - m*\} >5/2 
for ah {E*,m*) G M{Eo,mo,h)}. 

If e is sufficiently small depending only on 6 and the fixed parameters, then 
for any n, B{e,6,n) =^ F{e,n,As). 



Proof. Suppose that B{e,6,n) has happened. Recall that U{(j),e) is always 
non-empty by Lemma \20[ choose any U G ly({(p,e). Let E := Gh,n{4>iU'^) 
and m := Mh^n{4', U). Suppose that \E - E*\ < 6/2 and \m-m*\ < 6/2 for 
some {E* ,m*) £ M{EQ,mQ, h). We prove that this is impossible by arriving 
at a contradiction. 

Li the situation described above, the event K{e, n, E, m) has happened. 
Thus, there is a function r] satisfying the conditions (a), (b) and (c) of 
Lemma [22l Consequently, we have 

(al) ||(/)-ry||oo < e, 

(bl) 1(^0 - E*) - Hh,n{v)\ < 2e + 6/2, and 

(cl) |(mo - m*) - Mh,n{r])\ <Ce + Ce^'^ + 6/2. 
If e is sufficiently small depending only on 6 and the fixed parameters, then 
this shows that (/> is a (5-soliton, giving the desired contradiction. D 

Proof of Theorem [T71 Fix 6 E (0,1). Choose e E (0,1) small enough to 
satisfy the criterion of Lemma [25l Then for all n, 

F{B{e,6,n))<¥{F{e,n,As)). 

By the continuity of and the fact that As is a closed set not intersecting 
the region where G attains its maximum in TZ{Eo,mo, h), it follows that 

@{Eo,mo,h,As) < Q{Eo,mo,h). 
Lemma [Ml now completes the proof. D 
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12. Exponential decay of solitons 

Fix n. Suppose that q E C^" is a ground state soliton, that is, it minimizes 
energy among all functions with a given mass. Then by standard Euler- 
Lagrange theory, q satisfies 

-ujq= -Aq- \q\^~'^q, 

for some real number w. This can be rewritten as 

(67) {2dh-'^ + Lo)q{x) - h~'^ ^ q{y) = \q{xW~^ q{x) for all x € K- 

y ■■ y^x 
Let m := M/i,^(g). The following theorem shows that q must be exponen- 
tially decaying outside a small set. 

Theorem 26. There exists a subset U ofVn, whose size can be bounded by 
a number depending only on m, h, p and d, such that for all x G Vn, 

\q{x)\ < Ae-^^^(^\ 

where Djj{x) is the f^ distance of x from U, that is, the minimum of\y — x\i 
over all y G U, and A, b are positive constants depending only on m, h, 
p and d. Here y — x means the difference of y and x modulo n in each 
coordinate. 

The proof of Theorem [26] follows closely the outline of the proof of expo- 
nential decay in the continuum case, as given in [62\ Proposition B.7]. The 
main difference is that in the discrete case, we have to deal with discrete 
Green's functions. The proof is divided into several lemmas. 

Lemma 27. Ifp < 1+4/ d, then for any m > 0, h > and n > C{p, d, h, m), 

E^in{m, h) 

CO > ^ u. 

m 



Proof. Multiplying both sides of ([67|) by q{x) and summing over x £ Vn, we 
get 



to 



5^|g(x)|2 = -2d/i-2^|g(x)|2 + /i-2 ^ {q{x)qiy) + q{x)q{y)) 

x&V„ 



xeVn x£Vn ^,yev, 



xGV„ 



[x)\P+' 



= -h'' Y. k(^)-9(y)l'+El'?C 

x,y€Vn xeVn 



V " 

x.yeVn ^ ' X&Vn 



= -2h '^E^i^{m,h,n). 

By Lemma [9l Emini'm',h,n) — )• £'mm("T'5^) as n — )• oo and by Lemma [TT| 
EniiniiT^, h) < 0. This completes the proof. D 
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Lemma 28. Let r := 2d/{2d + u}h?). Letp{x,y,k) be the probability that a 
simple symmetric random, walk on the torus Vn starting at x at step is at 
y at step k. Then the soliton q satisfies for all x €z Vn the identity 

yev„ k=o 



Proof. Given q satisfying (I67p . let / be the function on the right-hand side 
in the above display. (Note that the series converges because r < 1, by 
Lemma [271) Our goal is to show that q = f ■ First, note that 

z : zr^x y£V„ fc=0 z : Zr^X ^ 

By the translation invariance of the torus, it is easy to see that 
^ p{z,y,k)= ^ p{x,w,k). 

z : z^x w : w^y 

Again, note that the random walk can be at y at time A; + 1 if and only if it 
was at some neighbor of y at time k and moved to y at the {k + l)th step. 
Therefore, 



p{x,y,k + l) = — ^ p{x,w,k). 

w : w^y 

Combining the last three displays, we get 



Y, fiz) = h^Y.Y.'^'^'p(^^y^^ + ^My)r\iy) 

z ■■ zr^x y&Vn k=0 

= y J2J2'-'^'p(^^y^''My)r\{y) 

y€V„ k=0 

-/i2^p(x,y,0)|g(y)rig(y) 

= {2d + ojh'^)f{x) - h'^\q{x)\P'^q{x). 
Comparing with (j67p . this shows that for all x G T^, 

{2d + uh^){fix)-q{x))- Y. ifiy)-q{y))=0. 

y : y^x 

In other words, (wl — A)(/ — g) = 0, where / is the identity matrix in 
C "^ ". Since A is a negative semidefinite operator (Lemma [7]) and u; > 
by Lemma [271 ojI — A is non-singular. This shows that q = f and completes 
the proof. D 

For any 5 > 0, let 

Us := {xeVn: \q{x)\ > 6}. 
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For X G Vn, let Ds{x) denote the £^ distance of x from Us, that is, the 
minimum of |x — y|i over all y £ Us- Here, as usual the difference x — y is 
computed modulo n in each coordinate. 

Lemma 29. For each x £ Vn and 6 > 0, 

\q{x)\ < Cor^'^''^+oj-^SP-^ ^ r^^^'^^qiy)], 

yeV„ 

where Cq = uj-'^mP+'^h-'^^P-'^'^d-^ and r = 2d/{2d + uh^). 

Proof. A random walk starting at x at time cannot reach y before time 
\y — x\i. Thus, p{x,y,k) = for all k < \y — x\i. By Lemma [28| this gives 

\Q{x)\<^^^\q{y)\p( Y^ r^+M 

y&Vn ^k=\x-y\i ^ 

yey„ 

Now, if y ^ Us, then \q{y)\P < dP~'^\q{y)\. On the other hand, if y e Us, 
then \y — x\i > Ds{x). But again, |g(x)p < h~ m for all x and 

\Us\ ^ un^. ^ 



Thus, 

Y r\''-y\^q{y)\P < -^r^'^'=\h-'^m)P + 6P~^ ^ ^'^"'''^(y)!- 
y&v„ y^Us 

This completes the proof. D 

Proof of Theorem [26l Define 

B{x) := Comaxr5l2^-^li+^^(?^), 
y&V„ 

where Cq is the constant from Lemma [29l Note that B{x) is never zero, 
B{x) > Cor^*(^' for all x, and for all x,y, 

B{y) = Comaxr5l^-S'li+^^(^) 

< Comaxrsl^'^li-ll^'^^li+^^W 
2ey„ 



\\y-^\^B 



r 2i« --i^i^ X . 



Let K be the smallest number such that \q{x)\ < KB{x) for all x G Vn- 
Since B is never zero on Vn and T^i is a finite set, K must be finite. By 
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Lemma [29| Lemma [271 and the above observations, 

\q{x)\ < Cor^'^""^ +Ku;-^6P~^ ^ rlJ'-^liS(y) 

yeVn 

< B{x) + KB{x)uj-^5P-^ Y^ r3ls/-^li 

yeVn 

< B{x) + KB{x)5^-^C{p, d, h, m). 

If 6 is chosen so small that 6^~^C{p, d, h, m) < 1/2, then the above inequality 
implies that 

- 2 

In other words, K <2. Thus, with such a choice of 6, 

\q{x)\ < 2Comaxr^l?'-^li+^^(2^) 
y£V„ 

for all X £ T^t- To complete the proof, note that for any y, 

-\y - x\i + Ds{y) > -\y - x\i H — 

^1, I , Ds{x) - \y- x\i Ds{x) 



2'" '2 2 

Thus, \q{x)\ < 2Cor^^(^)/2 for all j; G I4. D 

13. Lower bound 
Fix h > 0, and some Eq gM and mo > such that 

drriQ 



E^in{mo,h) < Eo < 2 • 

Let {E* ,m*) be a point in M{EQ,mQ,h). By Lemma I16[ m* is strictly 
positive. The numbers p, d, /i, Eq, mo, E* and ?n* will be fixed throughout 
this section and will be called the 'fixed parameters'. Any constant that 
depends only on the fixed parameters will be denoted simply by C, instead of 
C{p, d, h, EQ,m,o, E*,m,*). If the constant depends on additional parameters 
a,b,. . ., then it will be denoted by C{a,b, . . .). 

Recall the random function 4> defined in Section [6] and the objects M, TZ, 
Q and G defined in Section [lOl 

Theorem 30. Assume the condition (I68p . For arbitrary e € (0,1), let 
Bo = Bo{e,n) be the event 

{\Hh,n{^) - Eo\ < e, \Mh,n{^) - mo\ < e}. 

ThcTi 

liminfliminf'°«'^'^°'^-"»> 1 - ™, + §(£„,»„,/,). 

e^O n— s>oo n 
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Fix n and let /* be an element of C " such that Mh^n{f*) = "nio — m* 
and Hfi^nif*) = -E'mm("T'0 — m*,h,n). Then /* is a ground state soliton for 
the DNLS on Vn- Fix e > 0. Let C/ be a subset of Vn such that 

(69) M;,,„(r,C/^u5C/)<e2, 

where, as before, U^ = Vn\U. By Theorem [26l there exists a U satisfying 
the above property such that \U\ < C(e). 

Let z be a point chosen uniformly at random from Vn- Let (p' be an 
independent copy of (j) and define 

, , [^'(a;) if X e z + U, 
7(xj := < 

I (^{x) otherwise. 

Here z -\-U is the translate of U hy z on the torus Vn', that is, the addition 
is modulo n in each coordinate. 

Let U' be the set of all points that are either in U or adjacent to some 
point in U. In other words, U' = U L) dV^. Define the following events: 



Ai:={\Mh,n{(t>)-m*\<^}, 



A 



2 



Mh,n{(p) m* 



<e\ 



max|,/.(x)|2<n-'^/2^|0(x)|n, 

^3 := {MhA4>, z + U')< 4m*{2d + l)n-'^\U\}, 

Ai := {\(l)'{x) - f*{x - z)\ < n'^'^ for all x G z + U]. 

Finally, let A = A(e, n) := Ai n A2 fl ^3 PI A^. 

Lemma 31. Let A and 7 he defined as above. We claim that if A happens, 
then 

(70) \MhM - mo\ < 2e2 + C{e)n-'' 
and 

(71) \Hh,n{l)-Eo\<Ce^ + b{e,n), 

where h[e, n) is a number depending only on the fixed parameters, e and n 
such that lim„_>oo b{e, n) = 0. 

Proof. Suppose that A has happened. To prove (fTOj) . note that 

(72) Mh,n{7) = Mh.M, z + C/^) + Mh,n{cp', z + U). 
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By Ai and A3, 

\Mh,n{<t>,z + U'')-mr\ 

< \Mh,n{^) -m*\+ Mh,n{(p, z + U') 

(73) < e^ + C(e)n-'^. 
Now if A4 holds, then by the inequahty 

|a^-6n <r|a-6|max{a^-\6^-H < r\a - b\{a + \a - biy-^ 

that holds for any a,b > and r > 1, we see that for any x £ U and 
any r > 1, 

(74) ||r(x)r-|</.'(z + x)r| <C7(r)n-2'^. 

Thus by A4, (I69p . and the fact that Mh^nif*) = "t-o — "i*, we have 
|M;,,„(</.',z + [/)-(mo-m*)| 
< \Mh,n{r^ U) - (mo - m*)| + C\U\n-^'' 

(75) < e^ + C7(e)n-2«'. 
Combining ([72]), ([73D and ([75]) gives 

\Mh,n{l) - mo\ < 2e2 + C(e)n-'^. 
This proves (j70|) . Next, note that 

^h,n(7) = ^/^,n(</', ^ + t/') + i^h,n(0', ^ + U) 

xGz + U, y€z + U^ 

By dMD, dZlD and A4, 

5^ W{x)\''< Y, |r(x)|2 + C(6)n-2'^ 

< Ce^ + C(e)n-2'^. 
By^s, 

Y. |0(x)p < C(6)n-'^. 

The last two displays imply that 

(77) ^ Y. W{x)-ct>{y)\''<C^ + C{e)n-\ 
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By ^4 and ^, 

\Hh,n{(t>', z + U)- Hh,n{r^ U)\ < C(e)n-^\ 
Again from (j69p . it follows easily that 

From the last two displays, we have 

(78) |i^h,n(0', z^V)- Hh,n{n\ < Ce^ + C{e)n-^''. 

Next, note that 

By ^3, 

|G/,,„(0, z + U"")- Gh,nm < CMh,n{<P, z + U')< C(e)n-'^. 

Again by Ai and A2 and the fact that m* > 0, it follows that if e is sufficiently 
small (depending only on the fixed parameters), then 

\Gh,ni^)-E*\<Ce\ 

Lastly, note that by A2, 



Y^ \<p{x)r' < {mj,x\^ix)r') Y, \ 

xez + W X^Vn 



(X)|2 



The last four displays combine to give 

(79) \Hh,ni(p,z + U^) -E*\<Ce^ + C{e)n''^ + Cn-'^^P-^'^/^ . 

Combining ([76]), ([TTD, dTS]) and dZH) we get 

\HhAl) - (HhAn + E*)\< Ce^ + C{e)n-'' + 071-"'^^-'^^ 

By Lemma [9l H^^f*) — )■ £'min(?TT-o ~ "t-*) ^) as n — >• 00. On the other hand 
by Lemma [lU -Emin("io — f^* ■, h) = Eq — E* . This completes the proof. D 

Lemma 32. Let A = A{e, n) he the event defined immediately before the 
statement of Lemma [3T1 Then 

r r log^(-4(g,n)) ^ 

hm lim -j = 1 — 7710 + KyiEQ,mo, h). 

(Note that the definition of A(e, n) involves our choice of /*. The above 
result holds for any sequence of choices of /* as n — )• 00.) 

Proof. Let Ai, A2, A^ and A4 be as in Lemma [3ll Write 

P(^) = P(Ai) F(^2 1 ^1) P(^3 1 Ai n A2)F{A4 1 ^1 n ^2 n A3). 

For the first term, simply note that by Lemma [T8\ 

F{Ai) = exp(n'^ + n'^logm* - n'^m* +n'^o{e,n)), 
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where o(e, n) is a term depending only on e, n and the fixed parameters such 

that 

hmsuphmsup |o(e, n)| =0. 
e— >-0 rn-oo 

Next, define 






^1/2' 

Then, as in the proof of Lemma \2T\ S, is uniformly distributed on the unit 
sphere of C^", is independent of AIh,n{4')j satisfies 

and 

max.,,. |^(x)P ^ ^^ 

Consequently by Theorem 1151 and Proposition 1131 

-n'^^J — -^)+n'^o(e,n) 

Let z be as in the proof of Lemma [3TJ Note that since z is uniformly dis- 
tributed on Vn and is independent of (p, and therefore by Markov's inequality, 

E{Mh,n{<P,Z + U')\cl)) 



¥{Al I <P) < 



Am*{2d + l)n-d\U\ 



< 



^d Am*{2d + l)\U\ ~ 4m* 

Thus, if Ai happens and e is sufficiently small (depending only on the fixed 
parameters), then F^A^ \ cj)) > 1/2. Consequently, 

1 >P(yl3 I AinA2) > 1/2. 

Lastly, note that since the coordinates of (j)' are i.i.d. complex Gaussian with 
probability density function (nh) tt~^ exp(— (n/i) |xp), and are independent 
of z and (p, and \U\ < C(e); therefore by ([69]) . 

F{Ai I ^,z) =P(^4) = exp(-(n/i)'^^|/*(a;)|2+n'^o(e,n)J 

= exp(— n {rrio — m*) + n o{e, n)). 

Consequently, the same is true for P(A4 | Ai D A2 O A3). Combining the 
above estimates finishes the proof. D 

Proof of Theorem [30l By Lemma \3T\ if e is sufficiently small (depending 
only on the fixed parameters) and n sufficiently large (depending only on e 
and the fixed parameters), then the event A{e, n) implies the event BQ{e, n), 
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but with (j) replaced by 7. However, since 7 has the same distribution as (p, 
by Lemma [32] this completes the proof of Theorem [30l D 

14. The radiating case 

Fix /i > 0, and some Eq £M and mo > such that 

dmo 2dmQ 

As in Section [TTl the numbers p, d, h, Eq and rriQ will be fixed throughout 
this section and will be called the 'fixed parameters'. Any constant that 
depends only on the fixed parameters will be denoted simply by C, instead of 
C{p, d, h, Eo,mo). If the constant depends on additional parameters a,b, . . ., 
then it will be denoted by C{a, b, . . .). 

Recall the random function (p defined in Section [6j 

Theorem 33. Fix e G (0, 1) and 5 G (0, 1). Let Aq be the event 

{\Mh,n{(l)) - mo\ < e, \Hh,n{(l)) - ^o| < e}- 
Then for any fixed 5 £ (0, 1), 

logP(max^6y„ \(j){x)\ >5\Ao) 
hm sup hm sup ^ < 0. 

e— >0 n— >oo n 

The rest of this section is devoted to the proof of this theorem. Define 



?x:= 



H^) 



,1/2- 



{Ey^vMyry 

Then ^ is uniformly distributed on the unit sphere of C^" . Let 

2h'^Eo 



ao 



mo 



so that ao ^ [2d, id). Define three events: 

Ai := {\Mh,n{<P) - mo\ < e^}. 

A2:={ Yl le(^)-e(y)P- 



ao 






< e' 



yl3:={max|e(x)p<n-'^/2|_ 

xev„ 

Lemma 34. If e < C and n > C(e), then Air\A2C\ A'i 
Proof. Note that 



A.. 



E i^(^)-?(y)i'-«o 



x,y£Vn 
xr^y 



2h^ 



Gh,ni^) Eo 



Mh,n{<P) rno 



This shows that if e < C, then Ai D A2 implies that 
(80) \Gh,ni^) - Eo\ < e^/\ 
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Now, Aid A3 implies 

xeVn ^ 

X&V„ 



and hence 



xeVn 



Therefore if n > C(e), then A'^^,n(0) ^ e^'^. Combining this with 
completes the proof. D 

Lemma 35. 

lim inf lim inf -\ > 1 + log mo — mo — ^(^(ao)- 

e— >-0 n— >-oo n 

Proof. The random variable Mh^ni'P) ^'^'^ the random vector ^ are indepen- 
dent. Therefore by Lemma [Ml for e < C and n > C(e), 

(81) P(Ao) > F{Ai n ^2 n ylg) = P(Ai)P(yl2 n A3). 

By Theorem 1151 

(82) hm hm -j = -*rf(ao)- 

By Lemma [TSl 

[06) hm hm ;i = 1 + iogm-o — m-o- 

Combining (jSip . (j82p and ()83p proves the lemma. D 

Proof of Theorem [33l Define an event j44 as 

A4 := {max|(/)(x)| > 5}. 

x£V„ 

We have to show that 

^. ^. logP(AonA4)-logP(Ao) ^^ 
limsuplimsup ^ < U. 

If |-f^/i,n(</') — Eq\ < e and max^jgy^ |'^(a;)| ^ (J) then 



Gh,n{<t^) = Hh^M) + —-r V |0(x)| 
p + 1 ■^^ 
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Therefore, Aq n A^ implies 

|2 






X'~^y 



2h^{Eo-e + CdP+^) 



mo + e 
Thus, there is a constant Cq depending only on the fixed parameters such 
that if e < C{5), then Aq n A4 implies the event Ai n ^5, where 



■ x,yeVn 



Moreover, by the independence of Mh,n{4') ^"^^ ^i the events Ai and A^ are 
independent. Therefore, if e < C((5), then 

iP(^o n Ai) < P(Ai n A^) = P(Ai)P(A5). 

Theorem [TH shows that for any e < C{5), 

n— >oo n 

and we already know the limit of n~'^logP(yli) from (|83|) . Since ^d is 
a strictly increasing function in the interval [2d, Ad) by Proposition [T3l a 
combination of the above inequality and Lemma [351 proves the theorem. D 

15. Discrete concentration-compactness 

Fix h > and m > 0. Let /„ be a sequence of functions on Z such that 
Mhifn) = t^ for all n and Hh{fn) — >• £'min("T'j /i) as n — )• c«. The following 
theorem is the main result of this section. 

Theorem 36. There is a subsequence rik of natural numbers and a sequence 
of points yk € 'Z'^ such that the sequence of functions g^ defined as 

9k{x) ■■= fukiVk + x) 

converges to a limit function g in the L'^ norm for every q G [2, 00]. The 
function g has mass m and energy E^[n{m,h). 

The proof of Theorem [36] uses the well-known concentration-compactness 
argument (see [471 Section 1.4] and references therein). The main new chal- 
lenge in the discrete case is that the properties of E^i^{m,h) are not as 
well- understood as those of £"111111 ("i). 

For each x G Z and each positive integer R, let B{x,R) denote the i^ 
ball of radius R around x in Z . For any two positive integers n and R, 
define the 'concentration function' 



y£B{x,R) 



Pn{R) := sup }2 ^ l/-(^)l' 
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Note that Pn{R) is a non-decreasing function of R for each n. Let 

jj, := hm hniinf p„(i?). 

R-^oo n— !>oo 

Clearly, there are sequences R^ and n^ increasing to infinity such that 

lim PukiRk) = P- 

k—^oo 

It is easy to fix Rk such that Rk is even for each k. Passing to a subsequence 
if necessary (using a diagonal argument), we may assume that 

p{R) := lim pn,{R) 

fe— >oo 

exists for each positive integer R. 

Lemma 37. p = limfc^oo PukiRk) = hmfc^oo Pnk{Rk/'2) = lim/j^oo p{R)- 

Proof. First observe that from the monotonicity of /9„j. , 

limsup/)„j^(i?A:/2) < limsup pn^^iRk) = p. 

k—^oo fc— >oo 

On the other hand, for each R, 



and thus 



p{R) = liminf /9„ (i?) > liminf p„(i?). 



lim p{R) > p. 



Lastly, for any R, we have i?fe/2 > R for all k large enough, and thus 

PndRkm>Pn,{R)- 

Let k ^ oo gives that for all R, 

liminf p.„Ji?fc/2) > lim inf pn^{R) = p{R). 

Letting i? — )• oo completes the proof. D 

Since the function pn is bounded between and m for every n, therefore 
< p < m. 

Lemma 38. p > 0. 

Proof. Suppose that p = 0. Let all notation be as in Lemma [37l Then 
from Lemma Ell lim/j_!.oo p{R) = 0. But p is a non-decreasing, non-negative 
function. Therefore p{R) = for every R. In particular p{l) = and hence 
limfc__>oo Priki^) = 0. This implies that 

lim sup \fn^{x)\ =0, 

and therefore 



lim (Hhifn,) - Guifn,)) = lim — - ^ |/„,(x) 



x& 



<lim{snp\f^,ix)\Y~'^^^ = 0. 
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In particular, 

limmf Hhifn,) > 0, 

fe— >oo 

contradicting lim,fc_j.oo H^^fn^) = Eyain{fn^ h) < (Lemma [TT]) . D 

Lemma 39. There is a sequence of points yn in TL'^ such that 

liminf |/„(y„)| > 0. 

n— >oo 

Proof. Let all notation be as in Lemma [37l By Lemma [SU 

< /i = lim lim inf yO„ (i?) . 

Thus, there is some R such that 

liminf p„(i?) > 0. 

n— >oo 

Fix such an R. Let Xn{R) be a point such that 

PniR) = Y. h''\fniz)\'. 

Z&B{Xn{R),R) 

(The existence of such a point follows easily from the assumption that 
^hif^) = m < oo.) Since R is fixed, this shows that 

,. . fJ2z£B{x„{R),R)\fn{z)\ ^„ 
hmmt , , ' — > 0. 

\B{Xn{R),R)\ 

But there exists a point y„ G B{xn{R),R) such that 

., , ^|2 ^ ^z&B{x„{R),R) \fn(.z)\ 
\Jn\yn)\ ^ I „/ . jjs „N| 

\B{Xn{R),R)\ 

This completes the proof. D 

Lemma 40. If a and b are positive real numbers and a > 1, then 

(o + 6)">a° + 6". 

Moreover, for each < ci < C2, there is a positive constant C = C(ci, C2, a) 
such that whenever a,b £ [ci, C2], 

(a + b)" >a''+b''-C. 

Proof. The first inequality is a simple consequence of 

sup (x° + (1 - x)") = 1, 

0<x<l 

by putting x = a/{a + b). The second assertion follows similarly. D 

Lemma 41. For any positive m,m', 

Emm{m, h) + Erainim, h) > Erain{m + Ul' , h). 
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Proof. Let /„ and g„ be sequences of functions such that Mh{fn) = ni and 
Mh{gn) = fn' for all n, and 

lim Hhifn) = Emin{m,h), lim Hh{gn) = Emmim',h). 

Since for any function /, Mhi\f\) = M{f) and Hh{\f\) < Hh{f) by the 
triangle inequality, we may assume that the functions /„ and gn are non- 
negative real valued. By Lemma [39} there exist sequences of points y„ and 
Zn and a positive real number e such that 

liminf /„(y„) > e, liminf 5r„(z„) > e. 

n— >oo ra— >oo 

Define functions f„ as 

Vn{x) := /n(2; + Vn) + ^gn{x + Zn), 

where i = \/— 1. Since /„ and g„ are real valued, it is clear that Mh{vn) = 
m + m! and Gh{vn) = Ghifn) + Gh{gn) for each n. Now note that since 
p > 1, Lemma HOJ implies that for all x, 

bn(x)r+l = \fn{x + 2/n) + i5n(x + Z„)r+^ 

= ((/n(x + Vn)? + (5n(x + Zn)^)^^^^^^ 
>\Ux + ynW^^ + \gn{x + ZnW+\ 

If n is large enough, then fn{Vn) > e and gn{zn) > e- Since M/j(/„) = m 
and Mh{gn) = m' for each n, therefore the sequences fniUn) and gn{zn) are 
also uniformly bounded above. Therefore by Lemma 140} there is a positive 
constant C such that for all n, 

bn(0)r+l = ((/n(y„))' + (5„(Zn))')(^+')/' 
>l/n(yn)r' + bn(Zn)r'-C. 

Combining all of the above observations, it follows that 

\\Ta.suY>Hh{vn) < lim (Hhifn) + Hy,{gn))- 

n— >00 71— >CXD 

Since Mh{vn) = ra + m' for all n, this completes the proof of the lemma. D 

Lemma 42. n = m. 

Proof. By Lemma [38] we know that /x > 0. We also know from definition 
that < n < m. So we only have to eliminate the case < fi < m. Suppose 
that this is true. Let Xn{R) be as in the proof of Lemma [39l Then we can 
write 

fnk =Uk + Vk + Wk 

with 

Vk{x) = fnkix)l{\^-^^^(R^/2)\>R^}, 
Wk{x) = fnk{x)l{R^/2<\x-x„^{B^k/2)\<B^k}- 
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Then note that by Lemma [371 

xeBix„^ {Rk/2),Rk) xeBixr,^ {Rk/2),Rk/2) 

<Pn,{Rk)- Y. h''Wn,{x)\'' 

x&B{x„^{Rk/2),Rk/2) 

= Pnk{Rk) - PukiRkf^:) -> as A; -^ oo. 
By Lemma [37] it also foUows that 

Mh{uk) = Y h'^\fnkix)\'^ = PnkiRkf^) ^ P ask^(x,. 

x(^B{x„^^{Rk/2),Rk/2) 

Thus, linik^ao Mh{vk) = m — ji. From shnilar arguments using Lemma [371 
it follows that 

(84) lim {HhUn,) - {Hh{uk) + Hh{vk))) = 0. 

By the continuity of -Bmin (Lemma [9]) , 

liminfiJ/j(nfc) > Eminip,h), liminf i7/,(t!fc) > Eminim- fi,h). 

fc— >oo k—^oo 

Therefore by ()84p . 

liminf i:f,,(/„J > Eram{p, h) + Emin{m - fi, h). 

fc— >oo 

By the initial assumption that Hh[fn) — )• E^[^{m,h), this shows that 

Eminim, h) > -Emin(^, h) + E^i^{m - II, h), 
contradicting Lemma [JTJ D 

Proof of Theorem [36l By Lemma [321 we know that 

lim liminf yO„(i?) = m. 

_R— >oo n— s>oo 

Choose Rq so large that 

liminf p„(i?o) > rn/2. 
n— >oo 

Let Xn{R) be as in the proof of Lemma [39l Then for all sufficiently large n, 

Y h''\fn{x)\' > m/2. 

xeB{xn{Ro),Ro) 
Fix e G (0,m/2). Let R^ be so large that 

liminf /9n(i?e) > m — e. 

n— ^-oo 

Then for all n sufficiently large, 

Y /i'|/n(x)P>m-e. 

x£B{x„{Rc),Rc) 
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Since m — e + m/2 > m, and Mh{fn) = fn for all n, therefore for sufficiently 
large n, the balls B{xniRo), Ro) and B{xn{Re), Re) cannot be disjoint. In 
particular, 

\Xn{Ro) - Xn{Re)\l < Rq + Re, 

and therefore, 

B{Xn{Re),Re) Q BiXn{Ro) , Rq + 2Re) ■ 

Thus, if we define 

Vn{x) := fn{x + Xn{Ro)), 

then for every e G (0, m/2), there is a sufficiently large integer S^ such that 
for all sufficiently large n, 

Y, h''\vn{x)\^<e. 

X : |x|i>5e 

Thus, the sequence Vn is compact in L'^{'L ) and therefore has a convergence 
subsequence Vn^ , which we call gk ■ Let g denote the limit of gt in L^ . Since 
the convergence is in L^, it follows automatically that Mh{g) = m. Again, 
since the L°° norm of a function on Z'' is bounded above by its L^ norm, it 
follows that \\gk — g\\oo also goes to zero. Therefore, since for any q E (2, cxd), 

|2 



X] \Sk{x) -9{x)\'' < (sup \gk{x) - g{x)\y ^ ^ \gk{x) - g{ 
x&id ^^'^'^ xezd 



X 



it follows that gk converges to g in L*^ for any q £ (2,oo). This implies, in 
particular, that Hh{gk) — >• Hh{g). D 

16. Harmonic analysis on the lattice 

In this section p will not denote the nonlinearity parameter in the NLS. 
Instead, it will typically play the role of the p in the U* norm. 

We define the L'p norm for functions on U^ at grid size h as follows: 

It may seem strange to define a new norm by multiplying the standard 
l^p norm by a constant; the purpose of the definition is to ensure that the 
constants in the discrete analogs of classical inequalities (that we develop 
below) do not depend on the grid size h. Note also that ||/||p,/i = ||/||p, 
where / is the continuum image of / at grid size h. 

Similar inequalities were developed by Ladyzhenskaya |31j in the context 
of the "finite-difference method". However, since I could not find in |31j ex- 
actly what I needed (in particular, some delicate analyses of discrete Green's 
functions and a discrete version of the Hardy-Littlewood-Sobolev inequality 
of fractional integration), I decided to go ahead with my own derivations. 
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The grid size h will be fixed throughout this section. We will assume that 
h £ (0, 1), to avoid complications arising out of large values of h, in which 
we are not interested since we eventually want to send h to zero. 

16.1. Convolutions. We define the convolution of two functions / and g 
on Z at grid size h as 

Although the notation does not explicit include the grid size, it will be 
understood from the context. 

16.2. Young's inequality. Below we state and prove the discrete analog 
of Young's inequality for convolutions, at grid size h. The proof is exactly 
the same as in the continuous case; the important thing is that the constant 
does not depend on the grid size. 

Proposition 43. Let f,g be complex-valued functions on U^ . Let 1 < 

p^q^r < oo satisfy 

1 1 1 

- + - = - + 1. 
p q r 

Then for any h > 0, 

\\f *g\\r,h < \\f\\p,h\\9\\q,h- 

Proof Let a = (r — p)/r and f3 = {r — q)/r, so that a and j3 both belong to 

the interval [0, 1]. Let pi = p/a and p2 = <?//?, so that pi,P2 G [IjOo]. Note 

that 

1 1 1 

— + — + - = 1. 
Pi P2 r 

Let u = f * g. By Holder's inequality, 



\u{x)\ = h" 



J2 ■^(^ ~ y)9{y) 

Kh'^^lfix- y)|i-"|ff(y)|i-^|/(x - yr\giy)\'' 
<hHY.\f{x-yT'^>\g{y)\^^- ' 

\ i/pi / \ 1/P2 



yeZ'i yeZd- 

1/r 



h''Y.\f{x-y)\^'--y\g{y)\('~^>j ll/IISp„.ll5l,^,„,. 



yeZ'^ 

Taking rth power and summing over x gives 



1/9 



l7/ir < llflK^"")"" llnll^^"'^^'' llfll"'' llnll'^'' 
l«llr ^ \\J \\{l-a)r,hliyil(l-l3)r,hllJ Ilapi,h\\y''l3p2,h' 
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Since (1 — a)r = api = p and (1 — /3)r = I3p2 = q, this completes the 
proof. n 

16.3. Fourier transform. Take any h > and let K = 1/h. For a function 
/ : Z — 7- C with finite L^ norm, we define the Fourier transform of / at grid 
size h as 

where i = \/— T. Again, the notation does not explicitly include the grid size; 
it is to be understood from the context. It is easy to verify the inversion 
formula 

/(x) = 2-" [ f{i)e-'^^-< dC 

J[-K,KY 

With the above definition, \i u = f * g (at grid size /i), then 

provided that / and g are in L^ n L"^. Another easy fact is the Plancherel 
identity 

(85) II/II2,, = 2-^^ / Ifm'dC. 

J[-K,KY 

16.4. Littlewood-Paley decomposition. Fix h > and let K = l/h, as 
in the preceding Subsection. Let 70 : [—K, K] — )• [0, 1] be a smooth (i.e. C°°) 
function that is 1 in [—1,1] and outside (—2,2). Define 7 : [—K,K]'^ — )• 
[0, 1] as 

d 

i=l 

Note that 7 is simply g, where 51 : Z — )■ C is the function 
gix) = 2-'^ [ j{C)e-'^''^< dC 

J[-K,K]<i 
d ,-2 






where (/? : M — t- C is the function 

(86) (^(x) = \C 7o(t)e-'"* dt. 

For any a G (0,i^], let 7a : {-K.Kf -^ [0,1] be the function 7(^/0). A 
computation similar to the above shows that 7a is the Fourier transform of 
ga, where 

d 

ga{x) = a'^]^99(a7r/ixj). 
1=1 
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Lemma 44. For any a € (0, /f] and any p G [l,oo], 

Ci{p,d)a''^P~'^/P < \\ga\\p,h < C2(p,d)a'^(P-^)/P, 
where Ci{p,d) and C2{p,d) are constants depending only on p and d. 
Proof. Note that for any p € [1, oo], 



\9a\\ln = h''a''HY,W{a^hx)\ 






x& 



Prom the formula (j86p . the properties of 70, and standard results about 
oscillatory integrals (see e.g. [58, Chapter VIII, Proposition 1]) it follows 
that '^{x) decays faster than \x\~'^ for any a as \x\ — ?■ 00. Moreover, ip \s a 
continuous function. Therefore, as ah — t- 0, 

avr/i^ |(/7(a7r/ix)|P -^ / W{u)\^ du £ {Q, 00). 



x&L 

This shows that 



^l(P) ^ V^l ( J. MP ^ ^2(p) 

-^<\^\^i^a.hx)\^<^^.. 



which completes the proof. D 

Por each a, let Paf be the function whose Pourier transform is 

(7a(e)-7a/2(6)/(e)- 

2/2' ■ ' 



In other words, Paf = {ga — 9aj2) * /• Now, for any nonzero ^ G \—K^KY., 



00 



^{l2-iK'\i) - 72-0+1)^(6) = 1> 
i=0 

which shows that for functions with suitable decay at infinity (e.g. functions 

with bounded support) 

00 

j=0 

This is the discrete analog (at grid size /i) of the classical Littlewood-Paley 
decompositions (see e.g. [62l Appendix A]). 

Lemma 45. For any h > 0, any 1 < p < q, any a G {0,K], and any 
f :Z'^ ^C, 

\\Paf\kh<C{p,q,d)al-l\\f\\j,,h. 

Proof. By Young's inequality (Lemma |43]) . for any 1 < p < q and any 
aGiO,K], 

\\Paf\\q,h < (llSalU/i + ||9a/2 lU/i) ||/||p,/i, 
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where 

1 1 1 

- = - + 1--. 
s q p 

The proof is now easily completed using Lemma 031 D 

16.5. Gagliardo-Nirenberg inequality. The goal of this subsection is the 
prove a version of the Gagliardo-Nirenberg inequality for the lattice at grid 
size h. (For the well-known continuum version, see e.g. |62l Appendix A].) 
The result is stated as Proposition 38] below. To prepare for this, we first 
need some definitions and lemmas. 

For i = 1, . . . ,d, let €{ denote the ith coordinate vector. Let Vj denote 
the discrete derivative operator in the ith. coordinate direction, defined as 

(87) V./(x) := f(- + ^^l-f(-\ 

Note that Vj/ is simply the convolution of / with 6i, where 6i is the function 

'-h-'^-^ ifx = 0, 
6i{x) = < h^'^^^ if X = —Ci, 

for all other x. 

Lemma 46. For any ^ G [—K,K]'^, 

d 

Y.\vrm)\'>c\c\Vio\\ 
1=1 

where \S,\ is the Euclidean norm of ^ and C is a positive universal constant. 
Proof. A straightforward verification shows that 

v./(e) = — ^ — /(O- 

There is a positive constant Cq such that for all 9 G [— 7r,7r], 

|l-e^*''|>Co|0|. 
Therefore, if ^j G [-K,K], then 

\vrf{o\>cMfm. 

This completes the proof of the lemma. D 

Lemma 47. For any a G (0, K] and any f G L?{'L'^), 

\\Pafh,h < Ca-^Ghif), 
where C is a universal constant. 
Proof. If 1^1 < a/2, then |^j| < a/2 for each i, and hence 

7a(0 = 7(e/«) = 1 = 7(2C/a) = 7a/2(0- 
On the other hand, if |^| > 2vua, then |^j| > 2a for some i, and therefore 

7a (e) = 7(e/a) = = 7(2C/a) = 7a/2(0- 
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Combining these facts with Lemma H6t Lemma HS] and the Plancherel iden- 
tity (f85]l . we get 

\\Paf\\lh= I \K~f{i)?di 

J\-K.KY 

\{la{i)-la/2{i))f{0?di 

< ca-'- 1 E !(>(«) - 7a/2(e))vj(e)p d^ 

J a/2<\i\<2\/da ^^^ 
d d 

This completes the proof. D 

The fohowing proposition may be called the Gagliardo-Nirenberg inequal- 
ity for the lattice with grid size h. The important thing, as usual, is that 
the constant does not depend on h. 

Proposition 48. Take any 2 < q < oo and let 9 G (0, 1) solve 

I I 

Then for any f G Li{Z'^) D L'^{Z'^), we have 

\\f\U,h<C{q,d)\\f\\l^'G,{ff/'. 
Proof. Then for any /, by Lemma [45} 

oo 

\q,h < 2^ \\P2-3Kf\\q,h 
j=0 

oo 

<C{q,d)Y,{2-^K)-^-l\\P,-,Kfh,h 

3=0 

oo 

= C{q,d)Y,{2-^Kf\\P^-,Kf\\2,h. 
j=0 

Again by Lemma H5l 

(90) \\P2~.Kf\\2,h < Wfhh, 

and by Lemma H71 there is a universal constant Cq such that for all j > 0, 

(91) \\P2-^Kf\\2,h < CQ2lhGh{ff'\ 

Let Jo G Z be the unique integer such that 

Co2^o-i/iG;,(/)i/2 < ||/||2^^ < Co2^«/iG;,(/)i/2. 
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Then note that 

(92) V 2~'^K' < C{q, d)2-'^'K' < C{q, d) f ^'^f'^' ) 

and similarly 

J<Jo 

1-8 



Combining (l89]) . pID . (I9TD. (I92]l and (l93]l we get 



,h<C(g,d)5^(2-^K)lP2-.x/ll2,^ 



J=0 

00 

<C(g,d) ^ (2-^-i^)^||/||2,, 
j=io+i 
io 
+ C(g,d) ^ (2--'i^)^2^/jG,(/)V2 

j=-oo 

<C{q,d)\\f\\l^'G,{ff'. 
This completes the proof of the discrete Gagliardo-Nirenberg inequality. D 

16.6. Green's function. Let A be the Laplacian operator on C with grid 
size h, i.e. 

^/(^) = /|2 E (/(y)-/(^))- 

y ■ y^^ 

Lemma 49. Let I be the identity operator on C^ and co be any positive 
real number. For any u € L^(Z ), the unique solution to the equation 

{ul -A)f = u 

is given by f = g * u, where g is the discrete Green's function 

i2—d °° 

(94) g^^) = __Y,r''-'M^,k), 

fc=0 

where r = 2d/{2d + ujh'^) and p{x, k) is the probability that a d- dimensional 
simple symmetric random walk started at the origin is at x at time k. 

Proof. Note that A is negative semidefinite, since 



(/'^/) = -i E \f--fy\ 



■> I I — I 

2/i2 
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Thus for any positive w, wl — A is a positive definite operator. In particular, 
given a function u G L^(Z ), there can be at most one solution of 

(a;/ - A)/ = u. 

To show that g * u is a solution, one proceeds exactly as in the proof of 
Lemma m □ 



The Green's function is an indispensable tool in classical harmonic anal- 
ysis. While the continuum Green's functions are relatively simple objects, 
the discrete ones are more complicated. The purpose of this Subsection is 
to derive some careful estimates for the discrete Green's functions. 

Lemma 50. Let p{x, k) be as in Lemma\^[ Then for all x G Z and k > 0, 

p{x,k) < C(d)e-l^l'/2'=(l + fc)-'^/2^ 

where we interpret \x\'^ /2k as oo if x ^ and k = 0, and as if x = 
and k = 0. 

Proof. Suppose we are given x = [xi, . . . , Xd) £ Z*^, ^ > 0, and ki, . . . ,kd 
summing to k. If we know that for each i = 1, . . . , d, the walk has taken a 
total of ki steps along the ith coordinate axis, then the number of ways that 
the walk can be at x at step k is exactly 

n((fc,+L^)/2J' 

where the combinatorial term is understood to be zero if ki + Xi is odd, or 
if Xi [—ki,ki]. Therefore, 

p{x, k) = {2d)-^ y , , ^' n r ^\,] 

^^ ki\---kdli-i-\(ki + Xi) 2J 

0<fci,...,fc^<fc ^ "■ i=l ^^ ' 1-11/ 

fclH \-h^ — k 

fclH Vk^ — k 

Suppose k balls are dropped independently and uniformly at random into d 
boxes. Let Ki be the number of balls falling into the ith box. Then for any 
< ki, . . . ,kd <: k such that ki -\- ■ ■ ■ + k^ = k, 

kil---kdl 
Therefore, 

<'^' ''(-')=<n(,K,f;,)/2)^-")^ 
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A simple computation using Stirling's formula (e.g. the matching upper and 
lower bounds in [50]) shows that there is a universal constant C such that 
for any integers a > 1 and b G [—a, a] such that a + 6 is even, 



2 V a. 

It is easy to verify that for any x E (—1, oo), log(l + x) > x — ^. Applying 
this to control the logarithms on the right-hand side of the above expression, 
one gets 

\ g-fe2/2a 

]2-^<C- 



Sa + b)/2j - V« 
This holds for a > 1. To include a = 0, one modifies the inequality slightly 
to get 



(97) , ,,, ]2'''<C , 

^ ' \{a + b)/2) - TTT^ 

Using this bound in ([96]) and applying Holder's inequality, we get 
p{x, k) < C(d)E (e- Sti -V^K, ]^(i + i^,)- V2^ 

< C(d)e-l^l'/2^Em(l + Kir^'A 

(98) <C(d)e-l^'l'/2fcf J]E(l + i^,)"'^^M • 

(The interpretation for A; = is as in the statement of the Lemma.) 

Now, each Ki has a Binomial distribution with parameters k and Xjd. 
Therefore by Hoeffding's inequality [25] . 

P(i^i < kl2d) < e-^('^)^ 
which clearly shows that for any r > 0, 

E(l + Kiy < F{Ki < k/2d) + (1 + k/2d)-'' 
^^^^ <C{d,r){l + kr'^. 



Plugging this into (|98p proves the lemma. D 

Lemma 51. Suppose that f : [0,oo) — )• [0, oo) is a non-increasing function 
and g : [0, oo) — )• [0, oo) is a non- decreasing function. Then 

oo 



y^f{k + l)g{k)< / f{t)g{t)dt. 
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More generally, 

oo . 

Y.f{\x\ + Vd)g{\x\)< / f{\yM\y\)dy. 

Proof. For the first part, simply observe that when t E [fc,/c + 1], f{t) > 
f{k + 1) and g{t) > g{k), so that f{k + l)gik) < J^^' f{t)g{t)dt. 
For the second part, let Z+ = {0, 1, . . .} and M+ = [0, oo). Then 

oo 

Y, fWA + ^d)g{\x\) < 2^ ^ f{\x\ + ^d)g{\x\) 

and 

' f{\y\)g{\y\)dy = 2'' f fi\y\)gi\y\)dy, 



so that it suffices to prove 

oo „ 

^/(|x| + ^/dM|x|)< / f{\y\)g{\y\)dy. 



Take any x G Z^. Let B{x) be the set 

{y eM.'^ ■.Xi<yi<Xi + l, i = l,...,d}. 

The sets B(x) are pairwise disjoint and all have volume 1. Moreover, on the 
set B{x), \x\ < \y\ < \x\ + V", and therefore 

f{\x\ + Vd)g{\x\)< f f{\y\)g{\y\)dy. 

Jb(x) 

This completes the proof of the lemma. D 

Lemma 52. Given w > 0, let r = 2d/{2d + uoh?), as in Lemma H^ Take 
any /? > 0, a G M and x e Z'^. Let 

oo 
A:=0 

The following bounds hold: 

(a) If a > 2, then 

S < C{d,a,l3){Vd+ |^|)2-"e-c(d,/3,-)'^N. 

(b) If a = 2 then 

5 < C((i,a,/3)(C(/3) + C(/3)|log(/i(V^+ |x|))|)e-^('^''^'")^l"l. 

(c) If a < 2, then 

In all three cases, C{d,/3,uj) is an increasing function ofco. 
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Proof. First, assume that x / 0. For any y > 0, 



ry 1 
log(l + y)= -r^dz > 

Jo i + ^ 



1 + 2/ 
Consequently, for any K > and x £ [0, 1], 



l + isTx 
Since h E (0, 1), the above inequality shows that 

(100) r = — = i < e-^('^'^)'^' 

where C{d,uj) is an increasing function of oj. 

By Lemma [50] , Lemma [51] and the inequality (jlOOp , 

JO 

Applying the change-of- variable u = \x\'^/t in the above integration, we get 

where A = C{d,uj)h'^\x\'^. The inequality a? + h"^ > 2ab shows that for all 
u> 0, 

2n 2-2 
Consequently, 

S < |x|2-"e-^ TexpT-A _ ^')^{"-4)/2^^. 

When a > 2, the integrand may be bounded by a constant that depends 
only on /3, q and d (and not x), by simply dropping the A/2u term. 
Next, consider the case a = 2. Then 

Jya \ 2u 2 J J^ 

<C(/3) + C7(/3)|log^|, 

and by the change of variable z = A/u, 

[^ f A /3u\ _i , f°^ e-^/2 
/ exp \u du < I dz 

Jo ^v 2u 2 ; - J^ z 

<C + C\logA\. 
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Finally, consider the case a < 2. Again by the change of variable z = A/u, 



oo 

expl 



\ 2u 2 J - Jo 

< C7(a)^("-2)/2. 

This completes the proofs of all three cases when x ^ 0. When x = 0, the 
steps are essentially the same, but simpler. D 



Proposition 53. Suppose that d > 3. Given uj > 0, define g as in (I94p . 
Then for all x G Z , 

g{x) < C{d)h^-'^{Vd + \x\f-'^e-^^'^''^^^\''\, 

where C{d,u;) is an increasing function of uj (when d is fixed). 

Proof. First, suppose that x ^ 0. Then by Lemma H9] and Lemma [50} 

oo 

g{x) < C(d)/i2-'^^r'=+ie-l^l'/2^(l + A;)-'^/2. 

A:=0 

Lemma [52] now completes the proof. D 

Corollary 54. Suppose that d > 3 and g is the Green's function defined 
in (f94l) . Then for any t £ [l,d/{d — 2)), \\g\\t^h ^ C{t,d,uj), where C{t,d,uj) 
is a decreasing function of lo . 



Proof. Let A = 2\fd. By Proposition! 

ii9iiu<^' E c{tA){p'-\^rd^\x\f-'')' 



\x\<Alh 



E C'(t,d)e- 



Co(rf,w)/i|x|t 



\x\>Alh 



where Co(d, w) is an increasing function of lo. By Lemma [5T] and the as- 
sumption that /i < 1, 

\x\>Alh 



< ^Co{d,u^)tVd^d I e-^°^'^^^^''\y\'l{\y\yA/h}dy 



^Co{d,u;)tVd I " ^d-l^-Co{d,Lo)ut^^^ 
J A 



Since A = 2v«, it is easy to see that the last expression above can be 
bounded by a constant C(t, d,oj) that is a decreasing function of oj. 
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Again by Lemma ISTl 

\2-d\t 



h^ Y^ C{t,d){h^-'^{Vd + 



x\ 



'i / (h\„\\('^-d)t 



\x\<A/h 

= C{t,d) / u^-'+^^-^^'du. 

Jo 

Since h < 1 and d — 1 + (2 — d)t > — 1 + e for some positive e = e(t, d), this 
shows that the integrand is bounded by C{t, d) and completes the proof of 
the lemma. D 

16.7. Hardy-Littlewood-Sobolev inequality. Take any w > and let g 

be the Green's function defined in (j94p . Given a function u G L^(Z'^), let / 
be the unique solution to 

{ul - A)/ = u, 

so that by Lemma H9l f = g * u. The following proposition is a discrete 
analog of the Hardy-Littlewood-Sobolev theorem of fractional integration 
(see [57, Chapter V, Section 1.2]). 

Proposition 55. Let f and u be as above and suppose that d > 3. Let 
1 < p < q < oo satisfy 

1 1 2 

q p d' 

Then \\f\\q,h < C{p,q,d,uj)\\u\\p^h, where C{p,q,d,uj) is a decreasing function 
ofu. 

This theorem is actually quite a bit simpler than the classical Hardy- 
Littlewood-Sobolev theorem, which includes the endpoint case 1/q = 1/p — 
2/d. Including the endpoint requires a somewhat delicate argument using 
the Marcinkiewicz interpolation theorem, which we can afford to avoid. 

Proof. Let s = (0, 2) satisfy 

1 _ 1 s 
q p d' 

Let t = d/{d — s). Then by Young's inequality, 

\\f\\q,h < \\g\\t,h\\u\\p,h- 

Since t < d/{d — 2), therefore by Corollary [SU \\g\\t,h < C{p,q,d,u}), where 
C{p,q,d,uj) is a decreasing function of oj. This completes the proof. D 
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16.8. Derivatives of the Green's function. Recall the discrete derivative 
operator Vj defined in (187p . The following proposition gives an estimate on 
the size of Vj^f, where g is the discrete Green's function defined in (I94p . 

Proposition 56. Given cj > 0, define g as in (|94|) . If d > 2, then for all 
X G Z*^ and all 1 < i < d, 

where C{d,u;) is an increasing function of uj. When d= 1, 

|Vi5(x)| < C{1 + |log(/i(l + |x|))|)e-^('^)'^l^-|, 

where C{oj) is an increasing function ofuj. 

To prove this proposition, we first need to introduce some notation. For 
each non- negative integer k, let 

Pk := {{h,. ..,kd)eZ'^:0<h,...,kd<k, h + --- + kd = k}. 

For any k > and any (fci, . . . , kd) G Pk and {xi, . . . , x^) G Z , let 

'^^''■^''''■-''^-^'''--'''^■■=kj^.jl{i^ 

Here, as usual, we interpret (^) as if either a or 6 is not a non-negative 
integer, or if 6 > a. In other words, for the above expression to be non-zero, 
it is necessary and sufficient that for all j, Xj has the same parity as kj and 
satisfies \xj\ < kj. 

Lemma 57. For any k > 1, any 1 < i < d, any (ki, . . . ,kd) G Pk-i and 
any (xi, . . . , Xd) G Z'^, we have 

tp{k - l;ki,...,kd;xi,...,Xd) 

{ki + Xi + 2)d 



ip{k', ki, . . . , ki -\- I, . . . , kd'jXi, . . . ,Xi -\- i, . . . ,Xd)- 



k 



Proof. If kj + Xj is odd for some j, then both sides are zero. So assume that 
kj has the same parity as Xj for each j. Similarly, if \xj\ > kj for some j ^ i, 
then both sides are zero. So assume that \xj\ < kj for all j ^ i. 

If \xi\ > ki, then the left side is zero, and there are three possibilities for 
the right side. First, Xj may be equal to —k^ — 2. In that case, ki + Xi + 2 = 
and so the right side is zero. The other possibilities are that Xi < —ki — 2 
01 Xi > ki + 1. In both of these cases, |xj + 1| > /cj + 1, and hence the right 
side is zero. Note that the case Xi = —ki — 1 is excluded because Xi has the 
same parity as ki. Therefore in all three cases we have equality of the two 
sides. So we may now safely assume that |xj| < ki. 

At this point, we have that both sides are non- zero. Verifying the identity 
is now a simple algebraic exercise. D 
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Proof of Proposition [56l Take any x E Z such that '}2i=i^i i^ even. Then 
p{x,k) = for all odd k and p{x + ei,k) is zero for all even k, where e^ 



denotes the ith. coordinate vector. Thus, from the expression (j94p . 
g{x + ei) - 5((x) h^-'^ 



J2 r''^^p{x,k)- J2 r''^^p{x + ei,k)] 



k even fc odd 

(101) =^ E ^'""'(i-OK^;,^) 

A; even 



+ -^ Yl r''+\p{x, k-l)-p{x + e„ k)). 

k odd 

Now fix some odd k and some 1 < i < d. Then by the formula (j95|) . 

(102) p(x,A;-l)= ^ ■4){k-l]ki,...,kd;xi,...Xd) 

{ki,...,ka)ePk-i 

and 

(103) p(x + ej,/c)= ^ V(^;^i,--- ,fcd;2;i,...,Xi + 1,... ,xd). 

(fcl,...,fcd)GPfe 

Let Pj^ be the set of all {ki, . . . , kd) G Pk with ki / 0. It is easy to see that 
the map 

{ki,...,kd) ^ {ki,...,ki_i,ki + l,ki+i,...,kd) 
is a bijection between Pk-i and Pj^. Thus, by Lemma [571 

p(a;,A;-l)= ^ V(fc - 1; ^i) • • • > ^(i;a;i, • • • jXd) 

(A:i,...,A:d)ePfc_i 

= X] V'(^;^i,--->fcj + l>---,^d;2;i,---,a::i + !,•••, 2;^)^ ^^ 

{ki,...,ka)ePk-i 

= X] V'(^;fci,---,^d;a;i,---,a;i + l,---,a:d)^^ -^ • 

By (fT02]l and (fT03]l . this shows that 

\p{x,k - 1) -p{x + ei,k)\ 



< X ^(/c;A:i,...,/crf;xi,...,Xi + l,...,Xrf) 
(fci,...,A;d)ePs, 



{ki + Xi + l)d 

-^ l{fc,7^o} - 1 



INVARIANT MEASURES AND THE SOLITON RESOLUTION CONJECTURE 79 



Let Ki, . . . , Kd be the random variables defined in the proof of Lemma [5Ql 
Then by the above inequahty and (f97l) . we get 



\p{x,k - 1) -p{x + ei,k)\ 



<IE| 2"^ 



{Ki + Xi + l)d 



k 



l{i^.^o} - 1 



n 



K. 



\\{Kj+xj + i{j^,})r. 



< C((i)e-l^l'/2fcE 



{Ki + xi + l)d ^ 



k 



{K^T^O} 



11(1 + ^.)-^/' 



Now observe the following: 

• By dMD, for aU r > 0, E(l + Ki)''' < C{d,r){l + k)-\ By Holder's 



inequality, this gives 

{K, + x, + l)d 



E 



k 



l{i^,^o} - 1 



.7 = 1 ^ 



{Ki+Xi + l)d 



k 



< {l + k)-'^/^iE 

Since \xi\/k < {\x\/Vk)k-^/^ 
{Ki + Xi + l)d 



^{K,^o} - 1 



d+u l/{d+l) 



E 



k 



-1 



{K^T^O} 



1 



d 



,^.-V^.^.(E 



K,d 



d+U l/{d+l) 

d+U l/{d+l) 



-1 



{Kry^O} 



1 



Since Ki is a Binomial random variable with parameters k and 1/d, 
it follows by Hoeffding's tail bound l25] that ¥{Ki = 0) < e'^^'^'i^ . 
Moreover, K^ < k. Thus, 



E 



Kid_ 



Kid 



< e-<^(d)k + IE 



d+l^ l/(a!+l) 

d+K l/(d+l) 



k 



1 



Again by Hoeffding's bound, for any r > 0, 
Kid 



E 



/fc 



1 



<C(d,r)A; 



-r/2 



Combining all of the above, we see that for any odd k, 

(104) \p{x, k-l)-p{x + a, k)\ < Ci(d)A:-('='+i)/2g-C72(d)|x|Vfe 
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Using this estimate and Lemma [501 to bound the right-hand side in (jlOlj) . 
and applying Lemma [52} we get that for d>2, 

|Vi5(x)| < C7(d)/i3-'^(^/d+ |x|)2-'^e-^('='''^)'^l^l 

+ C(d)/ii^"'(Vd+ 1x1)1- V^('^''^)'^l^l, 

where C{d,uj) is an increasing function of uj. This proves Proposition [56] 
when d > 2. 

When d = 1, we use (fTIH|l . (fTUT]) . (IM]) and LemmaEHlto get 

oo 



A:=0 

Lemma [52] now completes the proof. D 

Corollary 58. For any d and any 1 < i < d, 

\\'^ig\\i,h < C{d,uj), 

where C{d,u;) is a decreasing function ofuj. 

Proof. Follows easily from Proposition l58] and Lemma IST] D 

17. Regularity of discrete solitons 

Suppose that 1 < p < 1 + A/d. As in Section [16l assume that h G (0, 1). 
Take any ?tt, > and let / G L?'{'L ) be a ground state soliton of mass m 
for the discrete NLS on Z under grid size h. That is, / minimizes energy 
among all functions of mass m. By Theorem 1361 it is easy to see that at least 
one such function exists. Also, by elementary Euler-Lagrange techniques, 
it follows that there is an a; > such that / satisfies the soliton equation 
{ojI — A)/ = 1/1^-1/, where / is the identity operator and A is the discrete 
Laplacian at grid size h as defined in (|16p . 

The purpose of this section is to prove regularity properties of discrete 
solitons. The main goal will be to prove that the smoothness bounds remain 
uniformly bounded as the grid size goes to zero. This is necessary for proving 
convergence to continuum solitons. The proof follows more or less the sketch 
of the proof of regularity for continuum solitons (see [62l Proposition B.7]), 
but using the discrete estimates from Section [16] 

In this section, C will denote any positive constant that depends only on 
p, d and m. In particular, C will not depend on h. Moreover, we impose the 
additional condition that C is uniformly bounded as m ranges over any given 
compact subinterval of (0,cxd). We will call this the uniform houndedness 
condition. If C depends on additional parameters a,b, . . ., then it is denoted 
as C{a, b, . . .). 

Lemma 59. u > 1/C. 
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Proof. Following exactly the same steps as in the proof of Lemma [271 we 
arrive at the inequality 

Emm{m,h) 
m 
The uniform boundedness condition holds due to Lemma [TTl D 

Lemma 60. G/,(/) < C. 

Proof. Taking q = p + 1 and 

n d{p-l) 



2(p + l)' 
the discrete Gagliardo-Nirenberg inequality (Proposition I48p implies that 

wfCU < c\\f\\^^,^'^^'^'^G,{fr^^'y' 

< cGhifY^^'^y^. 

(Note that, since h > 0, f £ L'^ implies that / E L'^. Also, it is easy to 
verify using the condition 1 < p < 1 + 4/d that 9 £ (0,1). Lastly, note 
that the uniform boundedness condition on C is clearly satisfied.) Now, 
Hh[f) = -E'min(m,/i) < by Lemma [TTl Consequently, by the previous 
display, 

Gh{f) = Hr,{f) + Nh{f) 

< c\\f\\;xi^ < cGh{ff^p-^y\ 

Since 1 < p < 1 + 4/(i, or equivalently, d{p — l)/4 G (0, 1), this shows that 
GhU) <G. U 

Lemma 61. For all q £ [2,oo], ||/||g,/i < C{q). 

Proof. The case g = 2 is already known from the assumption that Mh{f) = 
m. First, assume that q S (2, oo). If d < 2, then for any such q, we can find 
9 G (0, 1) satisfying ([88l) . and therefore by the discrete Gagliardo-Nirenberg 
inequality and Lemma [60l it follows that ||/||g,h < G{q) for all 2 < g < oo. 

Next, suppose that d > 3. By the discrete Gagliardo-Nirenberg inequality 
(Proposition [48]) and Lemma [60l we have that ||/||g,/i < G{q) for all 2 < g < 
2d/{d-2). 

Define a sequence qk as follows. Let go •= ^djid — 2). Note that the 
condition p < 1 + \jd implies that go > P + 1- For each A;, let g^+i satisfy 

1 p 2 



Qk+i Qk d 

unless the right hand side is nonpositive, in which let qk+i = oo. If g^ = oo 
for some k, then this definition implies that qj = oo for all j > k. 
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Note that if qh is finite for all k, then for each k we must have 
1 _p'' 2 
Qk qo d 



- = ^--{l+p + ...+p'' 1) 



k / A k—1 



i=0 



But as /c — >• cx), 



This shows that qk must become infinity at some finite k. 

Next, note that q^ is an increasing sequence. This is easily proved by 
induction as follows. Suppose that qk > qo = 2d/{d — 2). If qt+i = oo, there 
is nothing to prove. Otherwise, by the condition p < 1 + 4/(i, 

1 1 _ 2 p-l 

Ik Qk+i d qk 

^2 {p-l){d-2) 
- d 2d 

Take any k such that qk < oo. Suppose we have proved that for allq £ [2, g^), 

(105) \\f\U,h < C{q). 

Choose arbitrary g € [qk, qk+i)- Then l<qk/p<q<oo, and 

1 1 n 2 

_> >---;. 

q Qk+l Qk d 

Take q' G [2, gfc) so close to qk that 

1 p 2 
q q' d' 

Let /o := \f\P-^f. Since / = (w/- A)-Vo and to > 1/C by Lemma[5l it fol- 
lows by the discrete Hardy-Littlewood-Sobolev inequality (Proposition [55]) 
that 



\q < C'('7,^')ll/o||g'/p 

= Ciq,q')\\fr^,<Ciq,q'). 

However, q' can be chosen depending only on q, p and d. Since the sequence 
qk increases to infinity, this proves by induction that ||/||g,h < C{q) for all 
qe{2,oo). 

(An important thing to note is that we crucially used the fact that the 
constant in Proposition [55] is a decreasing function of lo, in conjunction 
with Lemma [59] to conclude that C{q) may be chosen to depend on lo only 
through m and not on the actual value of w, and that moreover, C{q) satisfies 
the uniform boundedness condition.) 
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Next, consider the case q = oo. Take any r > d/2. Let r' = r/{r — 1), so 
that r' < d/{d — 2). Then by Corohary 154^ Young's inequahty, Lemma [59] 
and what we have already proved above, 



|oo,/i — lb * /o||oo,h < ||5l|r',/i||/o||r-,h — lbllr-',/i||/||pr,/i < C- 

This completes the proof. D 

Lemma 62. For any q € [2, oo] and any i < i,j < d, 

\\V^f\\,,h<C{q) and ||V,V,/||q,;, < C(g). 

Proof. Since the ground state soliton minimizes energy, and the function |/| 
satisfies Mh{f) = Mh{\f\) and Hh{\f\) < Hh{f) by the triangle inequality, 
therefore we must have Hh{\f\) = Hh{f). Consequently, ||/(a;)| — |/(y)|| = 
\f{x) — f{y)\ for each neighboring pair of points {x,y), which shows that / 
must be of the form f{x) = afi{x) for some constant a G C with |a| = 1 
and some function /i : Z"^ — )• [0, oo). Therefore, without loss of generality 
we will assume in this proof that / is a non-negative function. 

By Young's inequality. Lemma [611 and the observation that Vj is a con- 
volution operator, we have that for any q £ [2,oo], 

\\V^f\U,h = \\{V,g)*f%,h 
^ ' <\\V^g\\l,h\\f%,h<Ciq)\\V^gh,h■ 

Next, note that by the inequality \a^ — bf\ < Taax{pa^~^ , plf~^}\a — b\ for 
non-negative a and b (which follows by the mean value theorem), and fact 
that ||/||oo,h ^ C from Lemma [611 we have that for all x, 

|V,r(x)|<C7|V,/(x)|. 
In particular, by Lemma [6l] this implies that for all q G [2,oo], 
(107) \\^if%,h < C\\Vif\\g,h < C{q). 

By the commutativity of convolution operators, 

ViV,-/ = (V,5) * (Vj/P). 
Thus, by Young's inequality and (jl07p . we see that for any q G [2,oo], 

||V,V,-/||g,, < ||Vi5||i,h||V,r||,,, 
^ ^ < C(g)||Vi5lli,h. 

Applying Corollary [58] (in conjunction with Lemma [59|) to ()106p and (llOSp 
completes the proof. The uniform boundedness condition follows from the 
monotonicity of the constant in Corollary [58l D 
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18. Continuum limit of discrete solitons 

In this section we establish that discrete sohtons converge to continuum 
sohtons as the grid size goes to zero, provided that the nonhnearity is mass- 
subcritical. 

Assume that 1 < p < 1 + A/d. Fix /i > 0, m > and let / be a ground 
state soliton of mass m for the DNLS system (I17p on Z at grid size h. 
Let Qx{m) be the unique ground state soliton of mass m for the continuum 
system ([3]). Let / be the continuum image of / at grid size h, as defined in 
Section M 

Theorem 63. For any q S [2,cxd], 

L'^iLQxirn)) < C{p,d,m,q,h), 

where C{p, d, m, q, h) satisfies, for any < tuq < rrii < oo and any fixed p, 
d and q, 

(109) lim sup C{p, d, m, q, h) = 0. 

The same bound also holds for \Hh{f) — E^i^{m)\ (with the modification that 
there is no q). 

In the following, C will denote any positive constant that depends only 
on p, d and m, satisfying the uniform boundedness condition defined in Sec- 
tion [T71 that is, C remains uniformly bounded as m varies over a compact 
subinterval of (0, oo), with p and d fixed. If C depends on additional pa- 
rameters a,b, . . ., then it will be written as C{a, 6, . . .). We will also adopt 
the convention that o(l) denotes any constant that depends only on p, d, m 
and h and satisfies (I109|) . 

Let -u; : K — >• R be the function 

fl-\t\ if|t|<l, 
wit) := < 
^^ [0 if|t|>l. 

Extend w to M.'^ as 

w(ti, . . . , trf) := w{ti)w{t2) ■ ■ ■ w{td). 
Define a function f^ : W^ — t- C as 

r(y):=^/(x)^(^). 

Lemma 64. The function f^ is absolutely continuous. If x ^ If and y = 
hx + ht for some t G (0, l)*^, then 

d 
f^{y)= Y. f{x + s)J[t^{l-t,f-^^. 
se{o,i}'' «=i 
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If dif^ denotes the partial derivative of f in the ith coordinate, then for x 
and y as above, 

se{a,\)<i i<i<d 

Proof. Since w is absolutely continuous with bounded support, it follows 
easily from the definition of f^ that f^ is an absolutely continuous function 
on M . Take x and y as in the statement of the lemma. Take any z G Z , 
and note that w{{y — hz)/h) is non-zero if and only if \yi/h — Zj| < 1 for 
i = 1, . . . , d. Since yi/h = Xi + ti for some ti G (0, 1) for each i, therefore 
w{{y — hz)/h) 7^ if and only if each Zi is either xi or Xi + 1. In other words, 

se{o,i}'* ^ ^ 

d 
= Y f{x + s)W^w{ti- Si). 

se{o,i}'* «=i 

An easy verification shows that when si S {0, 1}, w{ti — Si) = t^'(l — tiY~^^. 
This completes the proof of the first identity. For the second, note that by 
the previous display, 

^^f%y) = \^f%hx + ht) 
n oti 

Now, 

d il ifsi = l, 

-7^w{ti -Si) = \ 

oti 1-1 if Si = 0. 

This proves the second identity. D 

Lemma 65. The function f^ satisfies 

\Mhif) - M{f')\ < Ch, and \Hf,{f) - H{f')\ < Ch'/(P+'\ 

Proof. For x = {xi, . . . , Xd) € Z'^, let B{x) be the cube 

B{x) := {yeM.'^:xi< y^/h <Xi + l, i = l,...,d}. 

Another way to represent B{x) is {hx + ht : i E [0, 1) }. If y is a point in the 
interior of B{x), then Lemma [Ml shows that /'^(y) is a convex combination 
of {f{x + s) : s G {0, 1}'^}. Therefore, for any r G [2,oo), 

(110) \ny)r< max i/(x+s)r< J2 i/(^+^)r- 

^ ' se{o,i}'^ 
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Consequently, 

Summing over x gives 

/ \r{y)\^-dy=Y, I \ny)Ydy 

<2V^|/(x)r = 2^||/IU. 

x6Z<* 

By Lemma [m the bounds piOp and (jllip show that for all r G [2, oo], 

(112) \\nr<C{r). 

Next, note that 

ir(y)-/(x)|< max |/(x + s)-/(x)| 
se{o,i}<* 

d 

(113) < ^ Y,h\V,f{x + s)\. 

Recall that if y G B{x), then f{y) = f{x). Therefore by inequality ()113p 
and Lemma [60l 



/ (r(y) - hy)fdy =Y. I ifiy) - f{x)fdy 
< Ch'GhU) < Ch\ 



It is easy to see that for all r, ||/||j. = ||/||r-,/i- Therefore, for any r G [2,oo) 
by the L°° bounds from (J112p and Lemma [6T] and the above inequality. 



liirii.-ii/iir|<iir-/ii. 

<\\r - fwl^'^'wr - fwT <c{r)h'/^'. 

Using this bound, and again applying the L^ bounds from inequality ()112p 
and Lemma inn we get 

(114) \\\r\\:-\\fYr\<c{r)h'/'-. 

Next, let dif^ be the partial derivative of f^ in the ith coordinate. By 
Lemma [Ml 

Idifiy) - Vif{x)\ < max |Vi/(x + s) - Vi/(x)| 
se{o,i}'' 

d 

(115) < ^ j;/i|V,Vi/(x + s)|. 

se{o,i}'* i=i 
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Let dif be the function that is identically equal to 'Vif{x) in the interior of 
the box B{x), and arbitrarily defined on the boundaries. Then the above 
inequality, together with Lemma [62} implies that 

/ (a,r(y) - dJiy)fdy = Y. [ ^^^f'^y) " VJ(x))'<iy 

< Ch^ max ||VjVi/||^^ < Ch'^. 

Clearly for all r, ||9j/||r = ||Vj/||r,/i. Therefore by the above inequality and 
Lemma[62l ||5i/'^||2 and ||9i/||2 are both bounded by C. Consequently, again 
applying the previous display, 

\\\din\i-\\dj\\i\<c\\\dj'^h-\\djh\ 
<c\\d,r-djh 

<Ch. 
Using the above bound and (|114p . we get 

|M,(/)-M(r)| = |||/||2-||/^||2|<c/., 
and 

\Hn{f) - mm < \\\djg - iia,ri|2| + \\\fXtl - mCW 

<C/i2/(p+i). 

This completes the proof of the lemma. D 

Let Q be the unique, positive radially symmetric solution of ([7]). Let 
Q' '■— Qx(m') be the continuum ground state soliton of mass m' := M(f'^). 
Let Q : Z'^ ^ M be the function Q{x) := Q'{hx). 

Lemma 66. Recall the o(l) convention introduced immediately after the 
statement of Theorem I63[ Then 

\M{Q') - Mh{Q)\ < o{l) and \H{Q') - Hh{Q)\ < o{l). 

Proof. The function Q is a Schwartz function (see the remark following the 
proof of Proposition B.7 in [621 Appendix B]). In other words, Q is a C°° 
function, and all its derivatives decay faster than any polynomial. This 
provides ample regularity of Q, and together with the scaling relation ([8]), 
this completes the proof of this lemma. D 

Lemma 67. |F(/^) - H{Q')\ < o(l). 

Proof. Let a be a number such that aQ has the same mass as /. In other 
words, 



a 



MhiQ)' 
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By Lemma EH |M/,(/) - M(/'=)| < o(l), and by Lemma [Ml \M{Q') - 
M{Q)\ < o(l). But by definition of Q' , M{Q') = M(/^). Thus, 

(116) |a-l|<o(l). 

Now recall that: 

• Q is a Schwartz function (see the remark in the proof of Lemma [661) . 

• |M(/'^) — m\ < o(l) by Lemma [HSl 

• Q' is related to Q by the scaling relation ([8]). 

Combining the above, it follows easily that for all r G [2, oo) and 1 < i < d, 
||<5||r,ft ^ C{r) and ||Vj(5||r,/i < C{r). Together with (|116p . this implies that 

\H{aQ)-H{Q)\ <o(l). 

Since / is a discrete ground state soliton, therefore Hh{aQ) > Hh[f). Com- 
bining this with the previous display and Lemmas [65] and [66l we get 

H[Q')<H{n<Ht,{f) + o{l) 

< HhiaQ) + o{l) 

< Hh{Q) + o{l) < H{Q') + o{l). 

This completes the proof of the lemma. D 

Lemma 68. For any q G [2, oo] and any absolutely continuous v : M'* — t- C 
such that V G L^(M ) and \/v is uniformly bounded, we have 

Hull < r(d n\\\->,\\'^+^^ «^ «IIViill''+2^ «^ 

where \\Vv\\oo := maxi<j<d ||9jf ||oo- 

Proof. If 1 1 Vu I loo = 0, then since v £ L'^, therefore v must be zero almost 
everywhere. So let us assume that ||Vw||oo £ (0, oo). 

Take any xq G K'^- Let B be the ball of radius r around xq, where 

,.|Vi;||ooy 
Note that 

f^\v{x)\^dx ^ \\v\\l 
Yol{B) - Yol{B) ' 
which shows that there exists y £ B such that 

Since ||xo — y\\ < r, this shows that 

\v{xo)\ < \v{y)\ + ||xo - y||||Vv||oo 

< C{d)r-'^/^\\v\\2 + Cr\\Vv\\oo 
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Since this is true for every xq, the right-hand side is a bound for ||w||oo- To 
complete the proof, note that for any q G (2,oo), \\v\\q < ||f||^ Iblli- '-' 

Proof of Theorem [63l By Lemma [65l 

(117) \M{f'')-m\ = \m' -m\<o{l). 

Again by Lemma EH \H{f^) - H{Q')\ < o(l). But H{Q') = E^i^{m'), 
\m! — m\ < 0(1) and -Emin is a continuous function by ([6]); consequently, 

(118) \H{f')-E^Um)\<o{l). 

Together with Lemma [65t this proves that \Hh{f) — £'jnm("^)| < o(l). 

The bounds (J117p and (jllSp , together with the orbital stability of ground 
state solitons (see e.g. [47^ Proposition 3]) imply that 

^'(r,QAM)<o(i). 

By (jllSp and Lemma [621 ||9j/'^||oo < C for each i. Therefore by Lemma [68| 
L^r^Qx^m)) < 0(1) for each q£[2,oo]. D 

Corollary 69. For every fixed < ttiq < mi < 00, 

lim sup \Emin{m,h) - Emin{m)\ = 0. 

h^^O rno<m<mi 

Proof. For every m > and h > 0, let /m,/i be a discrete ground state 
soliton of mass m for the DNLS at grid size h. Then E^i^^m, h) = Hh{fm,h)- 
Theorem 1631 completes the proof. D 

19. Continuum limit of the variational problem 

Recall the objects 0, A^ an 7^ defined in Section [TOj and the function ^d 
defined in (I19p . The following theorem shows that as the grid size goes to 
zero, the set A^(£'o, m-Oi ^) converges to the single point (0, Eq — £'inin('n'o))- 

Theorem 70. Fix rriQ > and E^[^{mQ) < Eq < 00. For each h > 0, let 
{Efi,mh) be an element of Ai{EQ,mo,h). Then 

lim ruh = and lim E^ = En — Ejamiinr,). 

Proof. Recall the functions E^ (m, h) and E^ {m, h). If /i is sufficiently small, 
then £'0 < dmo/h'^ and by Lemma [69l Eq > Eyain{miQ,h). Therefore by 
Lemma [T6| when h is sufficiently small. 

Eh = E~^{mh, h) = Eq- £'min("io - rrih, h). 

We will prove that mh — )• as /i — )• by a subsequence argument. Then 
by the above identity and Lemma [69| the limit for Eh will be automatically 
established. Let ruh tend to a point m' G [0, tuq] along a sequence hi — )• 0. 
Then by Corollary [69l 

hm Eh^ = E' := Eq - Emm{mQ - m!) > 0. 
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We will show that m' = by the method of contradiction. Suppose that 
m' > 0. Fix m" G (0, m') and let 

E" ■.= Eo-E^in{mo-m"). 

Note that since -Emin is strictly decreasing (by Lemma SH and Lemma [TTIl . 
therefore 

(119) E" > E' >0. 
Also, defining 

E" := Eq - £'min(mo - m" , hi) = E^{m", hi), 
we see that by Corollary [691 

(120) lim E- = E". 

Fix 7 G (0, 1) and note that for all i, 
Q{Eh^,mh^,hi) -2log hi 

< log Eh,- [ log( ^^~^^^^^^' +27Vsin^(7rXj))dxi---(ixrf, 
^[0,1]'* \ ^h, ^ / 

and therefore (since m' > 0), 

limsup(G(£;/j^,m/j^,/ii) -2loghi) < logE' -C{'y), 



1,2 1 It, 7 



:= / log! 27 > sin [irxj) ]dxi ■ ■ ■ dxd- 



where 

r / '^ 

Ch) 

Since this is true for all 7 G (0, 1), we can take 7 — )• 1 in the above bound 
and get 

(121) limsup(e(£;,,^,m,,,,/ii)-21og/i,) <logE'-C(l). 

On the other hand, for any i, (|12Up gives 
e{E'/,m",h^)-2loghi 

= log E'/ - sup / log ( ^ „ ' ' + 27 y^ sin^ (vrx,- ) ) dxi • • • dxd 

0<'y<lJ[0,l]<i \ m" ^ J 

> log E'^ - [ log f ^ + 2 V sin2(^x,)') dxi • • • dx^ 

^logE" -C{1) asi^oo. 
Therefore, by (|121|) and (|119|) . this shows that for all sufficiently large i, 

(122) Q{E'i;,m",hi) > eiEh,,mh„hi). 
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But by the definition of £"", we know that {E^,m") G TZ{EQ,7no,hi). This 
contradicts the definition of {Eh^jTiih-) as a maximizer of Q(E,m,h) in 
niEo,mo,h). ' ' D 

20. Proofs of Theorems El [3] and H 

Proof of Theorem^ Lenima[T2] gives the formula for E^g,:^{m, h). LemmafTT] 
and Lemma [8] show that — oo < E^[^{m,h) < 0. The subadditive inequal- 
ity follows from Lemma HU Lastly, Corollary [69] shows that E^i^{m,h) — )• 
Eram{i^) as ^ — )• 0, and that the convergence is uniform over compact subsets 
of(0,oo). D 

Proof of Theorem [3l Take any sequence of function /^ on U^ such that 
Mhifk) -^ "i and Hh{fk) — >• £'min(W';^) as /c — ;■ oo. Let q^ be a con- 
stant such that M{akfk) = "t- Then a^ — t- 1 and therefore Hfi{akfk) also 
tends to £^mm("^, /i). Theorem [36] now guarantees the existence of a subse- 
quence oikjfkj converging to a limit / G S{ni, h) in the L'^ pseudometric for 
every q G [2,oo]. Now, /^ is uniformly L? bounded, and hence uniformly 
L"? bounded for every q G [2,oo]. Therefore, since a^. — )• 1, /fc. also tends 
to / in L'i. This proves compactness of S{m,h). To prove that the set is 
non-empty, simply note that by the definition of E^i^{m,h), there exists a 
sequence fk satisfying M/j(/fc) = m for all k and Hh{fi?j — )• £"111111 (m, h). D 

Proof of Theorem UJ This is a direct consequence of Theorem [63] D 

21. Proof of Theorem [5] 

To be notationally compatible with the theorems of Sectionsll H ll3landll4 ] 
we will write Eq instead of E and rriQ instead of m. As in those sections, the 
numbers p, d, h, Eq and ttiq will be fixed throughout this section and will be 
called the 'fixed parameters'. Any positive constant that depends only on the 
fixed parameters will be denoted simply by C, instead of C{p, d, h, E^jTUq). 
If the constant depends on additional parameters a,b,..., then it will be 
denoted by C{a, b, . . .). 

Fix n, and recall the definition of 5-soliton from Section [TT] Recall also 
the random function (p defined in Section [6] and the objects 0, Ai and TZ 
defined in Section [10] Lastly, recall that the set S{m, h) denotes the set of 
discrete ground state solitons of mass m at grid size h. 

Given a function / : V^ — )• C, let /"^ : Z — )• C be the extension of / to 



yd 



defined as 




if 2; G Vn, 
ifx 0K- 

Having defined f^, say that / is an "improved (5-soliton" if there exists 
g-.Z'^^C such that 

(a) 11/'' -g||oo < <5, and 
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(b) there exists (E*,m*) G Ai{EQ,mQ,h) such that 
\{Eo - E*) - Hh{g)\ < S and 
\imo - m*) - Mhig)\ < 6. 

Lemma 71. If E^i^{mQ,h) < Eq < drriQ/h?, then the set K in the state- 
ment of Theorem [5] can he alternatively described as 

K = K{h) := {m' G [0, rri\ : m' = tuq — m* 

for some {E* ,m*) G M{EQ,mQ,h)}. 

Proof. Suppose that m! G [0, mo] maximizes 

q{m) := log(mo -m,)-^A 

\ tuq — m 

Let m* = mo — m' and E* = Eq — E-^aini'm' ,h) = E^{m*,h). Note that 
E^{m*,h) is necessarily positive, for otherwise q{m') would be — oo, which 
is impossible since it maximizes q{m) in the interval [0, mo] and q{m) > — oo 
for m sufficiently close to mo- This shows that {E*,m*) G TZ{Eo,mo,h). 
Now take any {E,m) G M.{Eo,mo,h). Then by Lemma [T6l E = E^{m,h), 
and therefore 

n^F h\ ^ ,T, f'^h^{Eo-E^i^{mo-m,h)y 

B{E, m, h) = log m — Wrf 



m 
= q{mo — m) < q{m') = @{E* ,m* , h). 

Thus, {E*,m*) eM{Eo,mo,h). 

Next, suppose that we are given m' such that m' = mo — m* for some 
{E*,m*) G M{Eo,mo,h). By Lemma[l6l 

Q{E*,m*,h) = q{m'). 

Take any m G [0,m,o]- We have to show that q{m') > q{m). If q{m) = — oo, 
there is nothing to prove. If not, then 

< E'o — Eynm{m, h) = E^{mo — m,h). 

Let mi := ttiq — m and Ei := E^[mo — m,h). The above display proves 
that 

max{£;" {mi,h),0} <Ei= E+ (mi , h) , 
and hence {Ei,mi) G TZ{Eo,mo,h). Thus, 

q{m') = Q{E*,m*,h) > e{Ei,mi,h) = q{m). 
This completes the proof of the lemma. D 

Lemma 72. For any r/ > 0, there exists 5 = S{t]) > depending only on 
rj and the fixed parameters (and not on n), such that if f : Vn ^ C is an 
improved 6-soliton, then there exists v G S{m',h) for some m' G K, such 
thatL°°{f'',v) <T]. 



INVARIANT MEASURES AND THE SOLITON RESOLUTION CONJECTURE 93 

Proof. We will argue by contradiction. Suppose that the statement of the 
theorem is false. Then there exists an 77 > 0, such that for every positive 
integer k, there exist n^ and a function fk : Vn^. — )• C such that fk is an 
improved A;~^-soliton, but L°°(/|,t;) > rj for all v £ Um'eK'Si'm' , h). 

For each k, let gk be a function on Z satisfying the requirements (a) 
and (b) in the definition of improved (5-soliton (with 5 = k""^ and / = 
fk)- Let {El,mD be the corresponding element of M.(Eo,7no,h). We will 
show that /^ approaches an element of [Jm'eK<S{Tn', h) in L°° pseudometric 
through a subsequence, which will give us the necessary contradiction. 

Since A4{EQ,mQ, h) is a compact set (Lemma ll6p . we may assume without 
loss of generality that {E^^^m*^) approaches a limit {E* ,m*) S Ai{EQ, itiq, h) 
as A; — )• 00. Let E' := Eq — E* and m' := ttiq — m* . Then Mh{gk) — ^ 
m' and Hh{gk) — )• E' . But by Lemma [El E' = E^naifn' -.h). Therefore, 
by Theorem [31 gk approaches some g G S{m',h) in the L°° pseudometric 
through a subsequence. Since ||/^ — ^fcHoo — 5- as A; — ;■ 00, this shows that /^ 
also approaches g in the L°° pseudometric through the same subsequence. 
But by Lemma ITTl m' G K. This completes the argument. D 

Given a function / : V^ — )• C, let /,- denote a random translate of /, that 
is, 

fr{x):=f{x + T), 

where r is uniformly distributed on Vn and the addition on the right-hand 
side is addition modulo n in each coordinate. 

Lemma 73. // / is a 5-soliton, then, provided that n > C{5), 
P(/t- is an improved 2(5-soliton) > 1 — Cn'^ . 

Proof. Let g : Vn ^ C he a function satisfying the requirements (a) and 
(b) in the definition of J-soliton. Then clearly /,- is also a (5-soliton, with g^ 
serving the role of g. Let dVn denote the boundary of Vn in Z,'^. Since r is 
uniformly distributed on Vn, it is easy to see that 

\dV I 
E{Mh,n{9r,dVn)) = ^—fMhM < Cn-\ 

I ^n] 

Therefore by Markov's inequality, 

(124) nMhAdr^dVn) > n-V2) < Cn-^/^. 

Now, Mh{gr) = Mh^niOr) and Nhig^) = Nh^nigr)- Also, it is easy to verify 
that 

\Gh{gl) - Gh,n{9r)\ < CMh,n{9r,dVn). 

Thus, if Mh^n{gT,dVn) < n~^'^, then /,- is an improved ^'-soliton, where 
5' = S + Cn~'^/'^. By (fT2i|) . this completes the proof. D 

Proof of Theorem^ First, assume that 

(125) E,nm{'mo,h) < Eq < -£^max("lo,/i)- 
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By Lemma[l2l £'jnax("i-0) ^) = 2dmQ/h?'. Therefore, we are in the setting of 
Theorem [301 Let 

S ■■= Se,h,n{^0,'mo) 

= {ve C^" : \Mh,n{v) - mol < e, \Hh,n{v) - ^o| < e}, 

as defined in Section[2j Let / = fe,h,n be a random function chosen uniformly 
from S. Let 

A := {v £ S : V is not a (5-sohton}. 
Then by Theorem [T7] and Theorem [30l 

,. ^. logP(,^G^)-logP(0G5) 
iimsuplimsup ^ < u. 

By Lemma [6l this shows that 

iim sup lim sup ^ < ^• 

e— >0 n— >oo '^ 

In particular, 

(126) lim lim P(/ is a 5-soliton) = 1. 

But by Lemma [73l 

P(/t- is an improved 2(5-soliton) > (1 — Cn^ ) P(/ is a 5-soliton). 

However, /,- has the same distribution as /. Combined with (jl26p . this 
shows that 

lim lim P(/ is an improved 25-soliton) = 1. 

Since this is true for any 6 > 0, Lemma [72] completes the proof of Theorem[5| 
for the L°° pseudometric, under the condition (|125p . 

To prove the result for general q G (2,oo), simply observe that for any 
such q and any v ■.Z'^ ^ C, 

n??"! WiAii < ii?)ii'j~2||,,||2 

U^'J W^Wq ^ Irlloo 11^1125 

and that ||/||2 < C by definition. 

When £"0 > 2-E'max(w.o! ^)) Theorem [5] is a direct consequence of Theo- 
rem [33] and Lemma [U D 

22. Proof of Theorem [T] 

As in Section [211 we will write Eq instead of E and mo instead of m. The 
convention about the notation C will also be the same. 
Let K{h) be defined as in ([T23]l . Fix q £ (2, 00] and 5 > 0. 

Lemma 74. Whenever h < C{q,5), for all v € ^m'eK{h)S{''^' ■: ^) 
where v is the continuum image of v at grid size h. 
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Proof. Take a sequence hk decreasing to 0. For each k, let Vk be an element 
of Um/gx(/ifc)5(m',/ifc). Let m'^ := Mh^{vk) and m^ := mo - m'^. By defi- 
nition of ffc, m^ G K{hk)^ and by definition of K^hf.), there exists E"/; such 
that {Ek^rrik) G M.{Eo,mo,h). Therefore by Theorem [70l limfc_!.oo "^fc = 0. 
Consequently, Mh^{vk) — )■ ttt-q. Since we know that v^ is a discrete ground 
state soliton for every k, therefore by Theorem [H 

lim L'i{Vk,Qx(rno)) =0. 

A simple argument by contradiction now completes the proof. D 

Proof of Theorem^ Fix q G (2,oo] and 5 > 0. Take any h and any v G 
U^/g^(/j)5(m, /i). Note that 

L'^{fe,h,n: Q\(mo)) < -^^(A.M' v)+L'>{v, Q x(mo)) 
= L''{feAn, V)+L^V, Qximo))- 

Thus, if L'i{v, Qx{mo)) < '^/2 for all v G U.^/g;^(/,)5(m', /i), then 

lP(^^(AAn, Qxirr.0)) >S)<r{ inf inf Z«(/,,,,„, v) > 5/2). 

m'GK(h) vGo{m,h) 

Now, -Eminl^io) < -Eq < oo by assumption. Therefore, Corollary [69l and 
Lemma [T2] show that for all sufficiently small h, 

Therefore by Theorem [5] and Lemma [TU for all h < C{q,6), 
(128) limsuplimsupP(L«(A,ft,„, QA(mo)) > S) = 0. 

This completes the proof of the ffist part of Theorem [TJ For the second part, 
ffist fix g G (2, cx)] and 6 > 0. Fix some /c > 0. Choose hk so small that 
hk < 1/k and (jl28p holds with /i = /i^. Given /i^, choose e^ so small that 
ejt < 1/k and 

limsupP(Z''(/,^,,,^,„, gA(mo)) > '^) < 1/^- 

n— >oo 

Finally, given hk and e^, choose n^ so large that n^hk > k and 
This shows that 

hm P(Z«(4,^,,„,, QA(,no)) > 5) = 0. 

K— 5-00 

Since such a sequence {ek,hk,nk) exists for any 6 > (with g fixed), one 
can extract a sequence that works simultaneously for all (5 > by a diagonal 
argument. In particular, there exists a sequence (e^, hk,nk) such that for all 
<5>0, 

limP(L-(4,,,,„,, Qa(™o))>5)=0. 
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But again, for any t; : M — )■ C and any q G (2,oo], we have the in- 
equahty (jl27p . Since the L^ norm of f^^h,n is uniformly bounded (because 

||/e,/i,n||2 = Mh,nife,h,n) S [mo - e, mo + e]), this completes the proof of the 
second assertion of Theorem [H D 
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